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Abstract. The Lieb-Thirring inequalities give a bound on the negative eigen- 
values of a Schrodinger operator in terms of an L p norm of the potential. These 
are dual to bounds on the F'-norms of a system of orthonormal functions. 
Here we extend these bounds to analogous inequalities for perturbations of 
the Fermi sea of non-interacting particles, i.e., for perturbations of the contin- 
uous spectrum of the Laplacian by local potentials. 
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1. Introduction 

The Pauli exclusion principle for fermions in quantum mechanics has no clas- 
sical analogue. One of its primary effects is the increase in kinetic energy that 
accompanies an increase in the density of such particles. Intuitively, this increase 
should be quantifiable in a manner similar to that predicted by the semi-classical 
approximation to quantum mechanics, and it is the aim of this paper to show that 
this can, indeed, be achieved in the case of density perturbations of an ideal Fermi 
gas. 

We begin with some definitions [25l Chapters 3 and 4]. The state of a finite 
system of N fermions of q spin states each (q = 2 for electrons) is described by 
a density matrix T, which may or not be pure. Associated with a state is a one- 
body density matrix 7 (a reduction of T) which is an operator on L 2 (M. d , C q ). The 
essential properties of 7 are that < 7 < 1 as an operator and that Trj — N. It 
is a fact proved by Coleman in [5] (see also Thm. 3.2]) that any 7 with these 
properties arises from some state T, i.e., no other restrictions on 7 are required by 
quantum mechanics. 

The electron density is p~/(x) = Tro 7(2;, x), in which "/(x,x) is a q x q matrix. 
The kinetic energy of the N particle system depends only on 7 and is given by 
Tr(— A)7 in units where h = 2m — 1 and with A = V 2 denoting the Laplacian. 

The semi-classical approximation for the kinetic energy is 

Tr(-A) 7 «^ sc (d) / Pl (x) 1+2 / d dx , 

with the constant 



(1.1) K sc (d) := 



2 

d ( d{2n) d \ d 



d + 2 Vgl^- 1 



The Lieb-Thirring inequality [37] states that there is a constant < < 1 such 
that 

(1.2) Tr(-A) 7 > r d K sc {d) f p 1 {xj l + 2 ' d dx 



for any one-body density matrix < 7 < 1. The original value |27j for was 
0.185 but it has been improved since then to 0.672 jfi] . and it is current belief that 
it equals one (for d > 3). This subject continues to be actively studied (see for 
instance the recent works [TJ [T^l HB] and the reviews [TTJ HJ). Note that the 
inequality (|1.2[) does not require Tr 7 to be an integer; it need not even be finite. 

We can turn the matter around and, instead of specifying 7, think of specifying 
a density p(x) and asking for the minimum kinetic energy needed to achieve this 
particle density. The Lieb-Thirring inequality above gives a universal answer to 
this question in terms of the semi-classical approximation. Here, we are implicitly 
using the fact that for any given function p(x) > with J Rd p(x) dx = N there is 
a fermionic TV-particle density matrix whose one-body reduced density matrix 7 
satisfies Pl (x) = p[x), see [2U1 Thm. 1.2]. 

It is important for many applications that the right side of the inequality (|1.2[) 
is additive in position space. If we partition R 3 into disjoint subsets, the right side 
is just the sum of the corresponding local energies. While this does not hold for 
the left side, it nearly does. The bound shows that there is some truth to this 
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approximate additivity. This additivity, or locality, played an important role in a 
proof of the stability of matter [26] . 

While the inequality (|1.2[) was the object of principal interest in [33] , the actual 
proof of (|1.2[) went via the Legendre transform of (jl.2l) with respect to p. This is 
an inequality about the sum of the negative eigenvalues (Ej) of the Schrodinger 
operator — A + V for an arbitrary potential V, namely, 

(1.3) ^2\E j \=Tr(-A + V)-<f d L sc (d) [ V(xf+ d/2 dx, 

where X_ := max{— X, 0} > denotes the negative part of a number or a self- 
adjoint operator X, 

< L4 > = wpB*t • 

and ?d is a universal constant independent of V. In fact, the relation between the 
constants in (ll.2[) and (II ,3|) is given by 

(p L sc (dj) P ' (P' K sc (d)f = l and (f d f ' (rtf = 1 , 

where p = l + d/2 andp' = 1 + 2/d; see )26M25], The duality between p and V, and 
between kinetic energy and Schrodinger eigenvalue sums is one of the important 
inputs in density functional theory [20) . 

A question that is not only natural but of significance for condensed matter 
physics is the analogue of (|1.2[) when we start, not with the vacuum, but with a 
background of fermions with some prescribed constant density po > 0. How much 
kinetic energy does it then cost to make a local perturbation Sp(x) 1 ! This time dp 
can be negative, as long as p + Sp > everywhere. We would expect that the 
semi-classical expression will guide us here as well and, indeed, it does so, as we 
will show in this paper. 

The principal difficulty that has to be overcome is that inequality (|1.2p was 
obtained in [27] by first proving f| 1 . 3[) . a route that does not seem to be helpful 
now. The picture was changed by a paper of Rumin [35] in which inequality ()1 .2[) 
was obtained directly, without estimates on eigenvalues. (The constant obtained 
this way is not, however, as good as the 0.672 quoted above.) We are able to utilize 
some ideas in [35] to help solve our problem. 

The first thing is to formulate a mathematically precise statement of what it 
means to make a local perturbation of an ideal Fermi gas. One could think of 
putting N electrons in a large box of volume v, computing the change in kinetic 
energy, and then passing to the thermodynamic limit v — > oo with po = N/v 
fixed. For this, appropriate boundary conditions have to be imposed. To avoid this 
discussion we pose the problem for an infinite sea with specified chemical potential 
p > 0. The chemical potential of the ideal Fermi gas is 

A 4 = —j- ^sc(d) p Q < . 

It is often called the Fermi energy and can be interpreted, physically, as the kinetic 
energy needed to add one more particle to the Fermi sea. 

We then look at the operator —A — p, in L 2 (R d , C q ), which, in our context, plays 
the role of — A in inequality (|1.2|) . The energy observable of a particle is now defined 
to be —A — p, which is negative for states in the Fermi sea and positive for states 
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outside the sea. The energy to create either a particle outside the Fermi sea or a hole 
inside the sea is positive. The (grand-canonical) energy of the unperturbed Fermi 
sea is Tr(— A — /i)IT~ = — oo, where II - denotes the projection onto the negative 
spectral subspace of —A — fi. Clearly, Tl~(x. x) — po/qtc (we will often not write 
the identity matrix lc<? for simplicity). Our interest is in the formal difference in 
energy between the state described by a one-body density matrix < 7 < 1 and 
the state described by H~ , and this is non- negative since the minimum total energy 
(given fi) is the uniform, filled Fermi sea. Our main result is a lower bound for this 
difference in terms of the semi-classical expression in all dimensions d > 2, namely, 



for some universal r^ that does not depend on 7. The trace in this expression might 
not exist in the usual sense, that is, (—A — /i)(7 — II - ) might not be trace class. 
This situation will be dealt with more carefully in the sequel. 

The Legendre transform of the right side of the inequality (|1.5[) will give us an 
inequality for the change in energy of the Fermi sea when a one-body potential V 
is added to —A — p. The positive density analogue of (| 1 . 3[) in dimensions d > 2 is 



Of course, (— A — /Lt)_ = — (— A — /i)II~. The quantity above is not necessarily trace 
class but the trace can, nevertheless, be defined; see Definition 12.21 below. 

The inequalities (|1.5|) and (|1.6[) are only valid in dimensions d > 2. In dimension 
d = 1, a divergence related to the Peierls instability [29] appears, and a Lieb- 
Thirring inequality of the form of (|1.5I) or (11.61) cannot hold for p > 0. This will 
be discussed in detail in this paper. 

Our main inequalities (|1.5p and (jl.6|) were announced and discussed in [7 . In 
particular, the constants r 3 ~ 0.1279 and r-i ~ 0.04493 were given. We do not 
expect them to be optimal and it is a challenge to improve them. One interesting 
case in which the sharp constant in (ll.5[) can be found is that in which p^(x) is 
required to be zero for all x in some bounded domain f2. In Section 12.41 we prove 
that if the integral on the right side of (|1.5p is taken only over f2, then r<j = 1 in 
this case, and this is obviously optimal. 

Our method to prove the inequalities (|1.5j) and (|1.6|) is rather general and it can 
be used to treat other systems. As examples we will also discuss in this paper Lieb- 
Thirring inequalities in a box of size L>1 with periodic boundary conditions, the 
case of positive temperature, and systems with a periodic background. 

The paper is organized as follows. In the next section we introduce some math- 
ematical tools allowing us to give a rigorous meaning to the Lieb-Thirring in- 
equalities (|1.5p and (|1.6p . Our main task will be to correctly define the traces 
Tr(-A-/z)(7-n~) andTr ((-A- /i + !/)_- (-A -//)_) in such a way that flU]) 



(1.5) Tr (-A 




(1.6) 
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and (|1.6[) become dual to each other in the appropriate function spaces. In Sec- 
tion !2.3l we consider the case of a weak potential tV with t -C 1, and we compute the 
second-order term in t of the left side of (|1.6I) . This will clarify the fact that there 
cannot be simple Lieb-Thirring inequalities at positive density in dimension d = 1. 
The proofs of all these results are provided in Sections [3] and [4] In Section 12.41 we 
consider the case of a density matrix 7 which vanishes on a given domain f2 and 
we derive a lower bound on the relative kinetic energy which involves the sharp 
constant K sc (d). In Section[5]we prove Lieb-Thirring inequalities in a box with pe- 
riodic boundary conditions. This allows us to investigate the thermodynamic limit 
and to extend our results to positive temperature. Finally, in Section [6] we discuss 
the extension of our results to general background potentials, with an emphasis on 
periodic systems. 

Acknowledgment. Grants from the U.S. NSF PHY-1068285 (R.F.), PHY-0965859 
(E.L.), NSERC (R.S.) and from the ERC MNIQS-258023 (M.L.) are gratefully ac- 
knowledged. R.F. would like to thank Ari Laptev for a stimulating discussion. 

2. Statement of the main results 

In this section, we provide the necessary tools to give a clear mathematical mean- 
ing to the inequalities (|1.5[) and (|1.6|) which we have announced in the introduction, 
and we state our main results. 

We fix a positive number /1 > and denote by 

(2.1) ET := 1(-A < //) 

the spectral projection of the Laplacian associated with the interval (— 00,/x), de- 
scribing a free Fermi gas in its ground state with chemical potential /1. As recalled 
before, the ambient Hilbert space is L 2 (M. d , C q ) where q is the number of spin states 
per particle (which is 2 for unpolarized electrons but which will be taken arbitrary 
in this work). The gas, described by the projection II - , has the constant density 

a f q\S d ~ 1 \ d, 

(2-2) p = q(2n)- d dp= ^\^. 

J\ P \ 2 <n d(27r) d 

The kinetic energy per unit volume agrees with the semi-classical formula 

(2.3) q (2iry d f \p\ 2 dp = (d/2) L sc (d) = K sc (d) (p Q ) 1+ ^ 

J\p\ 2 <» 

where the semi-classical constants K sc (d) and L sc (d) are given by (|l.ip and (|1.4[) 
above. 

2.1. Lower bound on the variation of kinetic energy. We consider a fermionic 
state, with one-body density matrix < 7 < 1 acting on L 2 (R d , C 9 ), which we 
think of as a perturbation of the reference state LI - defined before in (|2.1j) . We are 
interested in proving a lower bound on the kinetic energy (including ji) of 7, counted 
relatively to that of II - , of the form of (|1.5[) . To make sense of this inequality, we 
use as main variable Q := 7 — II which satisfies the constraint 

(2.4) -IT < Q < 1 -rr := n+. 

Our goal is to prove a lower bound on Tr(— A — fi)Q. Our first task will be to give 
a clear meaning to this quantity, in a rather general sense. The constraint (|2.4p 
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can also be written (Q + IT") 2 < Q + IT - . Expanding (Q + IT") 2 shows that (f2~4|) 
is equivalent to 

(2.5) Q 2 < Q++ - Q— 

where we have introduced the notation Q TT := WQW for t, r' G {±}- In partic- 
ular, we have Q ++ > and Q < 0. Furthermore there is equality in (|2.5I) if and 
only if 7 — Q + II - is an orthogonal projection. 

For smooth-enough finite rank operators Q, the following computation is justi- 
fied: 



Tr(— A - n)Q = Tr (n + (-A - ^QTT + II" (-A - /i)QIT 

= Tr|-A- M |(Q++-Q--). 

As we have seen, we always have Q ++ — Q > 0, hence Tr(— A— fi)Q > (changing 
the density of particles inside or outside of the Fermi sea costs a positive energy 
once /z, the energy of the Fermi level, has been subtracted). We now use this fact 
to give a general meaning to Tr(— A — /i)Q, in the sense of quadratic forms. 

Definition 2.1 (Relative kinetic energy). Let Q be a bounded self-adjoint operator 
such that | — A — ^\ 1 ^ 2 Q ±± \ — A — /il 1 / 2 are trace-clas^. We define 



(2.6) 



Tr (-A - n) Q := Tr| - A - /*|3 (Q++ - Q~) | - A - n 



If Q is a bounded operator such that ±Q ±± > 0, then we extend the previous 
expression by letting 

Tr (—A - fi)Q := +oo 

whenever | - A - /i| 1/2 Q++| - A - p,] 1 / 2 or \ — A — fi\ 1/2 Q — | - A - fj,] 1 / 2 is not 
trace-class. 

Of course we have Tro(— A — fi)Q — Tr(— A — /i)Q (the usual trace) when (—A — 
fi)Q is trace-class. The previous definition of the relative kinetic energy is inspired 
by similar ideas used in the context of the Dirac equation 9 and of electrons in 
crystals 0]. Later on we will be interested in estimating the kinetic energy of 
operators of the form Qy = 1(— A + V < (i) — U for a given potential V. In 
general we do not expect such operators to be trace-class when \i > (or even 
compact, see Remark 12.41 below) . 



Remark 2.1. When 7 = Q + H is itself an orthogonal projection, 7 2 =7, we 
have equality in (I2.5[) and we obtain 

|-A-mI 1/2 T 



(2.7) Tr (- A - M ) Q = Tr | - A - ^ Q 2 \ - A - ^l 2 = 

where ©2 denotes the ideal of Hilbert-Schmidt operators on L 2 (R d ,C 9 ). 

We are now ready to state our rigorous version of (|1.5[) . 

Theorem 2.1 (Lieb-Thirring inequality, density version, d > 2). Assume that 
d > 2 and fi > 0. Let Q be a self-adjoint operator such that —IT" < Q < IT + and 



In the whole paper we use the notation Q for the two operators Q +Jr and Q , and the 
notation for Q +_ and + . 
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such that | — A — /i| 1/ ' 2 Q ±± | — A — p^ 2 are trace-class. Then Q is locally trace-class 
and the corresponding density satisfies 

(2.8) p Q € i 1+ 3(R d ) + L 2 (R d ). 

Moreover, there exists a positive constant K (d) < K sc (d) (depending only ond>2) 
such that 



(2.9) Tr (-A - p)Q 



> K(d) [ ({p Q + ()Q { x )) 1+ l-{p Q Y+l-l±^{p )- d p Q { x ) 
with pq the constant density of the Fermi gas, given by 



dx 



We recall that a locally trace-class self-adjoint operator A is such that Tr | Y^xl < 
oo for every bounded function \ of compact support. In this case, the associated 
density pA is the unique real-valued function in L\ QC {E. d ) satisfying Tr(yAx) — 
Ji^W 2 PA(x)dx. 

Note that since Q > H in Theorem 12.11 we have Pq(x) > —po for all x € M. d . 
The function 

(2.10) 5T; c (p) := (po + p) 1+ ~d - (po) 1+ ~d - ^(Pof* P 

is non- negative and convex for p > —po- Hence the integrand on the right side of 
(|2.9j) is always non- negative. The function ST^°(p) behaves like p 1+2 / d for large p, 
and like p 2 for small p. Moreover, it satisfies the scaling property 

(2.11) 5T^{p)^p d ' 2 5Tr{ P p- 2/d ) 
and one has 



]imST2 c (p) = p 1 - 



uniformly on M. + . In the limit p — > (which is the same as po —> by (|2.2p ). the 
inequality (I2.9[) reduces to the usual Lieb-Thirring inequality [371 HH1 US] 

(2.12) V0<7<1, Tr(-A)7 > K(d) [ p 7 (x) 1+ l dx. 

The best constant in this inequality is smaller than or equal to K sc (d), the semi- 
classical constant defined above in (|l.lj) . Hence, K(d) < K sc (d) must hold. From 
the scaling property (12.111) , we know that the best constant in (|2.9p is independent 
of p > 0. However, the best constant for p = in the Lieb-Thirring estimate (|2.12[) 
is not necessarily equal to the best constant for (|2.9I) . The recent estimates [6] for 
the Lieb-Thirring constant in (|2.12j) do not a priori give any information on the 
positive density analogue (|2.9|) . 

The proof of Theorem [53] is detailed later in Section 13721 It uses the convexity of 
(5TJ C , to estimate separately the densities corresponding to the two diagonal terms 
Q and the two off-diagonal terms Q ±T . The estimate on the diagonal terms 
<2 ±=t is based on a new method which has recently been introduced by Rumin |35) . 
This estimate works similarly in dimension d — 1. The off-diagonal terms Q ±T are 
studied by a direct and explicit method which does not cover the case d = 1. 

There cannot be an inequality like (|2.9[) in dimension d = 1 for p, > 0. This 
surprising fact is due to a special divergence of the off-diagonal terms Q ±T at 
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the Fermi points (see Section l2~3l below for details). However, we can prove the 
following: 



Theorem 2.2 (Licb-Thirring inequality, density version, d = 1). Assume that 
d = 1 and fi > 0. Let Q be a self-adjoint operator such that H < Q < IL + and 

such that | — A — /i] 1 2 (5 ±=t | — A — pj 1 ^ 2 are trace- class. Then Q is locally trace- class 
and the corresponding densities satisfy 

(2.13) p Q±±e L 3 (R)nL 2 ( R ), f |fc| ; ^ WT |^(fc)| a ctt<oo > 

(VM+ifci)iog ^ 2 ^:| fc || j 

(where^ denotes the Fourier transform). Moreover, there exist two positive con- 
stants K(l) < K sc (l) and K'(l) > such that 

(2.14) 

Tr (-A - n)Q 

^ K ( 1 ) / ({Po + Pq++(x) + Pq--(x)) 3 - (p ) 3 - 3(p ) 2 (pq++(x) + Pq—(x))j dx 

i K (^+|fe|)log(^^) 
wii/i po the constant density of the Fermi gas, given by (|2.2[) . 

Note the logarithmic divergence of the function in the denominator, at \k\ = Z^/fl. 
Hence the last term is not bounded from below by L \pq+- +Pq-+ | 2 • In Section [2~^l 
below, we will see that, up to the value of the prefactors K(\) and ^'(1), this bound 
is optimal. In particular, the right side of (|2.14l) cannot be replaced by a constant 
times J r o~T^ c (pq). In the limit \i — > 0, the inequality (|2 . 14[) nevertheless reduces 
to the one-dimensional Lieb-Thirring inequality (I2.12[) . 

Remark 2.2. Let cp 6 i 2 (R d ,C g ) be any normalized, smooth enough function. 
Applying (|2.9|) or f|2.14f) to Q± = ± |n ± (p)(n ± (p| and using a simple convexity 
argument, we obtain the following Sobolev-like inequality: 

(2.15) / |p 2 ~ M | |^(p)| 2 >i^(d) / <57; sc (M 2 ) 
for all tp with J„ d \ip\ 2 < 1, and in any dimension d > 1. 

2.2. Variation of energy in presence of an external potential. In this sec- 
tion we study the dual version of our Lieb-Thirring inequalities (I2.9P and (|2.14p . 
expressed in terms of an external potential V (the variable dual to p) . We will give 
a rigorous meaning to (|1.6[) . 

Let V be a real-valued function satisfying 

(2.16) V € L 2 (R d ) n L 1+ 2 QR d ) ford>2 



(2.17) V e L 3/2 (R) + L 2 {R) 
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Under our assumption (|2.16|) . the operator —A + V is self-adjoint on H 2 (M. d ), by 
the Rellich-Kato Theorem [33]. In dimension d = 1, our assumption (|2.17p allows to 
define the Friedrichs self-adjoint realization of —A + V, by the KLMN theorem |33) . 
We now define 



(2.18) 



Qv-=Uy-n^ where U v := 1(— A + V < ft), 



as well as U y := 1 — Ii v . 

Remark 2.3. The real number fi could a priori be an eigenvalue of — A + V. Then, 
Theorems 12.31 and 12.41 below hold exactly the same if Hy is replaced by Hy + S, 
where S is an orthogonal projection whose range is contained in ker(— A + V — fi). 
In dimension d > 3, it is indeed known [15] that, under our assumption (12.16[) on 
V, the self-adjoint operator —A + V has no positive eigenvalue, thus \x is not is 
the point spectrum of —A + V. However, \i could be an eigenvalue of —A + V in 
dimensions d = 1 and d — 2. 



Similarly as in Definition 12. 1[ we can define a relative total energy as follows. 

Definition 2.2 (Relative total energy). Let R be a bounded self-adjoint operator 
such that | - A - fi + V\ 1/2 Ii v RIi v \ - A - fi + V\ 1/2 are trace-class. We define 
(2.19) 



Try (- A — /x + V)R := Tr | - A — /i + V| 2 (n+iffl+ - U v RUy) \ - A - fi + V\ 



If R is a bounded operator such that ±H V RU V > 0, then we extend the previous 
expression by letting 

Try (-A - fi + V) R := +oo 

whenever \ - A- i^i + V^^UyRUyl- A- ^l 1 / 2 c 
is not trace-class. 



-A-^n-iffl-i-A- 



1/2 



Since Qy is the difference of the two orthogonal projections Ily and II , we have 
at the same time 



-IT < Q v < n+ and 



n y < -Qv < n+. 



Hence both Tro(— A — fi)Qv and Try(— A — (i+V)Qv make sense by Definitions [2T] 
and 12.21 With our definitions we have 



(2.20) 

and 

(2.21) 



Tr (-A-n)Q^ 



Try(-A - fx + V)Q\ 



A - /i 



1/2 



Qv | - A-n + V\ 



1/2 



In the theorem below, we show that, under suitable assumptions on V, the two 
quantities (|2.20|) and (|2.21|) arc finite and that 



(2.22) 



Tr y (-A - fi + V)Qv = Tr (-A - fi)Q v + [ Vp Q , 



as expected. We also derive an estimate on Try (—A — fi + V)Qv which is the dual 
version of (1231) for d > 2. 
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Theorem 2.3 (Lieb-Thirring inequality, potential version, d > 2). Assume that 
fi > and d > 2. Let V be a real-valued function in L 2 (R d ) n L 1+d / 2 (R d ). 

• Both Qv | — A — /ip/ 2 and Qv \ — A — /i + V] 1 ^ 2 are Hilbert- Schmidt operators, 
hence (I2.20|) and (|2.21|) are finite. 

• The relative total energy Try (—A — fi+ V)Qv can be expressed as 
(2.23) Tr v (-A - fj, + V)Q v 



= min Tr (-A-^)Q + / V(x)pQ{x)dx 

-n-<q<n+ V JR d 

|_A- Al | 1 /=Q±±|_ A - Al | 1/2 ee 1 

The minimum in this formula is attained for Q = Qv ■ In particular, (|2.22[) holds 
true. 

• We have the inequality 
(2.24) Thy(-A- p + V)Q v 

f / 1+- , , d 2 + fi d \ 

> -L(d) \(V(x) - /i)_ +2 - M 1+ 2 H ^— ^2 y(x)J 

'» = ^((STP®) ,2 « 



and where K(d) is the optimal constant in 

We recall that the semi-classical constant L sc (d) is defined above in (|1.4|) . 
Let us comment on our result. We can formally write 

(2.25) Try(-A-/i + y)Q v " = " -Trf(-A-/i + y)_ + (-A- A t)_)- / 9 / V 

where po is the constant density of the translation-invariant state IT - , recalled in 
(|2.2j) . The first term of the right side is the formal difference between the total 
(grand-canonical) energy of the Fermi gas in the presence of the local perturba- 
tion V, and its total (grand-canonical) energy in the translation-invariant setting 
without any potential. The term p J Rd V, which makes sense under the additional 
assumption that V G L 1 (R d ), is also the first order term obtained by perturbation 
theory when the first term is formally expanded in powers of V. 
The semi-classical approximation of the right side of (|2.25[) is 

f ( i+- , , d 9 + d d \ 
L sc (d) \(V(x) - M)_ 2 - M 1+2 H t^— V(x)j dx 

and, up to the value of the multiplicative constant L(d), it is precisely the right 
side of our estimate (|2.24l) . Our result therefore says that the variation of energy 
obtained by including the potential V in the system is O(l) in the thermodynamic 
limit, and (j2.24j) provides a precise estimate in terms of the size of V. Since the 
term po J R d V is obtained via first-order perturbation theory, the semi-classical term 
on the right side of (|2.24[) is therefore an estimate on the validity of the first order 
approximation. 
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In Scction[52J we will render the formal equality (|2.25p more rigorous, by means 
of a thermodynamic limit argument. More precisely, we show in Theorem 15.31 that 



(2.26) Tr y (-A- p + V)Q v 
= lim (- Tr L 2 (CL) (-A L - M + V1 Cl )- + Tr i2(Ci) (-A L - - p V 

where — is the Laplacian on a box Cl = [— L/2,L/2) d with periodic boundary 
conditions. This will also justify our definition of the total free energy. A tool to 
prove (|2 . 26[) is to derive a Lieb-Thirring inequality similar to f|2 . 24[) . for a system 
living in a box with periodic boundary conditions (Theorem 15. 2p . 

The estimate (|2.24[) follows from the density estimate (|2.9|) and the variational 
principle (|2.23|) . by noting that 

Tr (-A-/x)Q+ / Vp Q >K(d)f ST° c ( Pq ) + f V p Q 

JR d JR d JR d 

for all — LT~ < Q < II + , by Theorem [2Tj Optimizing the right side with respect to 
Pq (keeping in mind that pq is pointwise bounded from below by —po), yields (|2.24l) . 
Similarly, if we assume that (j2.24j) is known, we can derive (|2.9p by choosing 

Hence (|2.24[) and (|2.9[) are dual to each other. 

In dimension d = 1, using the weaker lower bound (|2 . 14[) on Tro(— A — p)Q, we 
can prove the following result. 

Theorem 2.4 (Lieb-Thirring inequality, potential version, d = 1). Assume that 
p > and d = 1. Let V G L 3 / 2 (R) + L 2 (R) be a real-valued function such that 

Then all the conclusions of Theorem \2.3\ remain true, except that (|2.24p must be 
replaced by 

(2.28) Tr v (-A - p + V)Q V > -1,(1) J ((V(x) - p)l - /J + | pk V(x) \ dx 



with 



2 / 1 \ 1/2 
L (!) = o olTTTT ^ M 1 ) <md £'(1) 



3 V 3-^(1) / ~ 4£T'(1) 



We will see in Section 1231 below that it is not possible to take L'(l) = 0. 
When p — > the inequalities (|2 . 24|) and (|2 . 28[> reduce again to the usual Lieb- 
Thirring inequality [26l [27] which is the dual version of (|2.12p : 

f i+- 

(2.29) < Tr(— A + V)_ < L(d) / V(x)_ 2 dx. 
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Remark 2.4. In general Qv is not a compact operator. Indeed, it was shown by 
Pushnitski [30] (see also [31]) that the essential spectrum of Qv is 

\\s(n)-it i-5( M )-iir 



2 2 



where S(fi) is the scattering matrix associated to the pair (—A,— A + V). Hence 
Qv is not compact, unless S(fi) = 1. Similarly, one does not expect, in general, 
that (—A — /j,)Qv and (—A — /x+ V^)Qy are trace-class, rendering Definitions 12.1 
and 12.21 necessary. 



Remark 2.5 (Relation with the spectral shift function). The spectral shift function 
Cv(A) formally satisfies [3D] 

(2.30) f ( v (\)d\ = -Tr v (-&- V + V)Q V - po f V. 

J -oo JR d 

If V is in L 1 (R d ) and satisfies the assumptions (|2. 16|) or (|2.17l) . it is possible to 
define Qv as the (distributional) derivative of the right side with respect to /i. 

2.3. Second-order perturbation theory and the ID case. In this section we 
compute the variation of energy when a potential tV is inserted in the system, to 
second-order in t. In particular we will show that in the one-dimensional case d = 1, 
the constant L'(l) in the lower bound (|2.28[) cannot be taken equal to 0. 

The following result, whose proof is sketched in Section I4T31 is well known in the 
physics literature [29] . 

Theorem 2.5 (Second-order perturbation theory). Assume d > 1. Let V be a 

real-valued function in L 1 (R d ) n L°°(]R' i ) and fi > 0. Then, using the notation of 
the previous section, 

where 



-d 



dp 



IpI 2 <i \p-k\ 2 -\p\ 2 
| P -fc| 2 >i 



(2.32) 



47r|fc| "VIS 
|5 d - 2 | Z" 1 2X/T - 72 + |A:| 



log " v / i ' -^r r rf " 2 dr i/ d > 2 



2|fc|(2^) rf 7 to ^v/T^ 2 "- \k\ : 

In particular, when d — 1 f/ie constant L'(l) appearing in (|2.28|) mitst satisfy 
L'(l) > g/(127r). 

Our proof is valid under much weaker assumptions on the potential V, but 
we have not tried to optimize this. The divergence at |fc| = 2yTfl of ^i(-/^/JI) in 
dimension d = 1 is well-known, and it is sometimes called the Peierls instability [291 
Sec. 4.3]. When the interactions among the particles are turned on, the system 
becomes unstable because of the large number of possible electron-hole excitations 
between the two points ±2^/j2. A macroscopic deformation of the system can 
sometimes lead to the opening of a gap at the Fermi points [29l HH EH EZ] EH [38] . 
In higher dimensions, the second-order response function VE^ is bounded (this also 
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follows from our bound (|2.24|1 ). but it is seen to have an infinite derivative at 
\k\ = 2y/JI, a fact sometimes referred to as Migdal-Kohn anomaly |28[ 116] . 

We note that the semi-classical approximation to the left side of (|2.31l) satisfies 

1+- , i d 2 4- d d \ 

(2.33) lim — * 2 - L = - d ^±2l^ V 2 . 

v ' t->o t 2 8 P 

This proves that for d = 1, it is not possible to take £'(1) = in (|2.28p . since 
the response function diverges at the Fermi points k — ±2^//Z whereas the semi- 
classical second-order term stays finite. A closer inspection of the constants reveals 
that L'{\) > q/(12w) must hold, as stated in Theorem 12.51 
It is possible to calculate ^2 and ^3 exactly: 



(2.34) * 2 (|fc|) 




(2.35) * 3 (|*|) = 
Furthermore, we have the following recursion relation 

\qd-3\ rl r d-3 f \U \ 

(2 .36) ( *, fo, t > 3. 

which implies that is strictly decreasing for all d > 3 (whereas for d = 2, ^2 is 
constant on [0, 2] and strictly decreasing on [2, 00)). We deduce that 

\ad-2\ i-l d-2 d(d-\-9) 

Observe that in dimensions d > 2, perturbation theory predicts the same value for 
the constant L(d) as semi-classics does. This is not so surprising since the largest 
constant is obtained if V is supported close to 0, hence V is very spread out in 
x-space, which puts us in the semi-classical regime. 



Remark 2.6. As is detailed in Section 14731 below, the second-order perturbation 
of the energy arises from the first-order term in the expansion of Qtv- This term 
is purely off-diagonal (the corresponding {Qtv) ±:t vanish to first order in t). This 
emphasizes the fact that the absence of a Lieb-Thirring inequality in ID is due to a 
possible divergence of the off-diagonal densities Pq±t in Fourier space at |fc| = 

The corresponding first-order density is proportional to $ d* V. For potentials V 
whose Fourier transform does not vanish at the Fermi surface, this density decays 
slowly in cc-space, due to the lack of regularity of at \k\ = 

2.4. A sharp inequality. We state and prove in this section a lower bound on the 
relative kinetic energy needed to banish all the particles in a domain, f2, from the 
Fermi gas. This inequality involves the sharp constant K sc (d) and it is the positive 
density analogue of a result due to Li and Yau |19) . 

Theorem 2.6 (A sharp estimate for the energy shift). Assume d > 1. Let ft be 

an open subset of R d of finite measure. Let /1 > and denote, as before, II := 
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1(— A < For any fermionic density matrix such that 

0< 7 <l B ,\ n , i.e., 0< (/,7/W«) < / \f(x)\ 2 dx, V/eI 2 (l d ), 

JR d \n 

and such that \ — A — p\ 1 ^ 2 Q ±± \ — A — p] 1 ^ 2 are trace-class, with Q := 7 — II - , we 
have 

(2.37) Tr (-A- M )Q > ^ sc (d)Po +5 |fi|. 

The constant in this inequality is best possible. In dimension d > 2, applying 
Theorem 12. II and using that p 7 (x) = on fi, we get 

Tr (-A -p)Q>K(d) pl +d/2 \n\ 

+ K{d) ( (p^ 2 p^ 2 -^ P d l\p, ~ Po)) dx . 
Jm. d \n K as 

Here K(d) is not optimal but, on the other hand, the bound also quantifies the 
fact that p 7 cannot be equal to po close to the boundary because of the Dirichlet 
conditions. 

Proof of Theorem \2.6[ We have 

q— = n-(7 - n-)ir < n-i R<xn n _ - n~ = -n-i n n-. 

Using that Q ++ > 0, we get 

Tr (-A - p)Q = Tr I - A - M | 1/2 (Q ++ - Q~)\ — A - p\ x l 2 
>Tr(-A-/z)L /2 :M-A-/z)y 2 

= (2^)- d |fi| / (|p| 2 - M )_dp. 
Recalling the definition of K sc (d), we obtain the claim. □ 

3. Kinetic energy inequalities: Proof of Theorems 12.11 and 12.21 

3.1. Preliminaries. In this section we state and prove some preliminary results 
that will be useful in the proof of our main theorems. 

Throughout the paper we denote by JC = ©oo (resp. B) the algebra of compact 
(resp. bounded) operators on L 2 (R d , C 9 ). The usual norm of bounded operators 
is simply denoted by || • ||. We also denote by G p (for 1 < p < 00) the ideal of 
compact operators A on L 2 (M. d , C 9 ) such that Tr \A\ P < 00, endowed with its norm 
\\A\\ &p = (Tr 

In order to simplify the statements below, we introduce the following Banach 
space 

(3.1) X := |q = Q* G B : Q\ - A - p\ x ' 2 G 6 2 , 

l-A-pl^Q^l-A-pl^ee^, 
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endowed with its natural norm 



(3.2) \Q\ X :=\Q\+ Q\-A-p\ 1 ' 2 + | - A - M | 1/2 Q ++ | - A - 



I' 



1/2 



Si 



-a-mI^q-i-a-mI 



1/2 



Si 

For the sake of simplicity, we do not emphasize the dependence in /x in our notation. 
The space X has a natural weak topology which is the intersection of the ones 
associated with the spaces appearing in the definition ()3.2[) of ||-|^. Here Q n — Q 
in X means Q n Q weakly-* in B, Q n \ — A — fi] 1 ^ 2 — v Q\ — A — /J.] 1 / 2 weakly in 
6 2 and l-A-^^Q^l-A-fi] 1 / 2 | — A — At| 1 ' /2 Q ±± | — A — /i.| 1/2 weakly-* in 
&i. The unit ball of X is weakly compact for this topology, by the Banach-Alaoglu 
theorem. The following convex subset of X will play an important role: 

(3.3) IC:={QeX : -IT < Q < II+}. 

Our first result deals with the continuity of the map Q E X h-> pq in L\ oc {K l ). 

Lemma 3.1 (Operators in X are locally trace-class). We assume that p > and 
d > 1. Let Q be a self-adjoint bounded operator in X. Then, for every bounded 
function r\ of compact support, there exists a constant C n such that 

(3-4) \¥2v\gi < C v Mx ■ 

Hence Q is locally trace class and pq is well-defined in Lj L oc (]R ci ). 

Furthermore, the map Q £ X t— > r\Qr\ G 6i is weakly continuous: If we have a 
sequence {Q n } such that Q n — v Q weakly in X, then r\Q n r\ — ¥ r\Qr\ strongly in &±. 
In particular, pq n — > pq strongly in L 1 1 oc (IR £i ). 

Proof. We consider the spectral projection III := l( — A < max(l,2/i)), which 
localizes in a ball containing strictly the Fermi surface, and we denote by II2 = 1 — TEi 
its complement. Then we write Q = J^ fe l=1 „ TlkQTli and estimate each term 
separately. We start with H2QH2 which we treat as follows 

V^Q^n = v - I - A - ^Q++\ - A - p\i ^~rV 

|-A-m| 2 \-A-p\2 

where we have used that II2 = Il2n + . Since 77 and II2I — A — p]^ 1 ^ 2 are bounded, 
it is clear that the previous operator is trace-class. Furthermore, we know that if 
T n — 1 T weakly-* in 61 and K is compact, then KT n K — > KTK strongly in ©1. 
Hence the weak continuity follows from the fact that 77II2I — A — p]^ 1 ^ 2 is compact. 
For niQri2, we write similarly 

jTfiQI^ = r?n x Q\ - A - n\i ^~rv 

\-A-p\2 

and use that rylli € 6 2 , Q\ - A - p\ l > 2 € 6 2 and n 2 | - A - p\-^ 2 i] e K. The 
argument is then similar as before. Finally, for iliQIli, we simply use that 77II1 G 
62 and that Q is bounded. The rest follows. □ 

Remark 3.1. The previous proof does not use the fact that | — A — n\ 1 / 2 Q~~~ \ 
A — pi] 1 ! 2 is trace-class. 

The following says that finite rank operators are dense in X in the appropriate 
sense. 
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Lemma 3.2 (Density of finite rank operators). For every Q £E X , there exists a 
sequence Q n € X of finite rank operators, such that (— A.)Q n 6 B and 

Qn ~> Q strongly (that is, Q n f — > Qf strongly in L 2 (R d , C q ) for every fixed 
f e L 2 {R d ,C q )); 

(Q n -Q)\-A-v\ 1/2 



• lim 

n— >-oo 

• lim 

n— >oo 



= 0; 

A- Ai | 1 / 2 (g„-g) ±± |-A- M | 1 / 2 



= 0; 

61 



PQ n -> PQ strongly in L\ QC {R d ). 

Furthermore, if Q belongs to the convex set K defined in (|3.3I) . then Q n can be 
chosen in K, for all n. 

Note that operators Q € X are not all compact, hence in general \\Q n — Q\\ -ft 0. 

Proof. We start by approximating Q by a sequence of Hilbert-Schmidt operators 
Q n , with (Q n ) ±± € 61. Let us define the orthogonal projection P n := l(l/n < 
I — A — fx\ < n), which localizes in momentum space away from the Fermi surface 
and from infinity. We now define Q n := P n QP n . It is easy to verify that Q n is 
a Hilbert-Schmidt operator by choice of P n and, similarly that ((5 n ) ±=t are trace- 
class. We have P n — > 1 strongly in L 2 (R d ). Since Q is bounded, we obtain that 
Qn Q strongly. Also, it is well-known that when A £ 6 p for some 1 < p < 00, 
then P n AP n — > A strongly in & p . In particular, we have that P n Q\ — A — /i] 1 / 2 /*, = 
Q n \ — A — /il 1 / 2 — > Q\ — A — fi] 1 / 2 strongly in 62, using that P n commutes with 
I — A — fi] 1 ' 2 . The convergence of the trace-class terms is similar, and the strong 
convergence of pQ n in Lj oc (lR d ) follows from Lemma [3.11 Finally, we note that, 
since P n commutes with II - , Q n belongs to K. for all n, whenever Q is itself in K.. 

For a proof that Q n can itself be approximated by smooth finite rank operators 
in /C, see [TQl Theorem 6]. □ 

3.2. Proof of Theorem I2.lt kinetic Lieb-Thirring inequality for d > 2. In 

this section we prove Theorem 12.11 for /1 > (the case fj, = is well-known [26j 
H3HS]). Replacing Q by U U QU* where (U u f)(x) := p d ^f( y /JIx), it is easy to 
verify that (|2.9[) follows from the case fi = 1, which we will assume throughout the 
proof. Also we assume for simplicity that the number of spin states is q = 1 but 
the proof for the general case is identical. Finally, since the semi-classical energy 
difference 6Ti c (defined in (|2.10p ) is non- negative, the right side of our Lieb-Thirring 
inequality (|2.9|) is lower semi-continuous with respect to pq. This shows that, by 
Lemma 13. 2\ we can prove (|2.9p assuming that Q is a smooth-enough finite rank 
operator, and deduce the general case by density. 

Recall our notation Q = n~QII~, Q ++ = n + QII + and so on. We will 
estimate the density arising from each term separately. The constraint — II~ < 
Q < n+ is equivalent to Q 2 < Q ++ - Q — . 



Step 1. Estimate on Q . In order to bound the density arising from the diagonal 
terms, we will use the following generalization of the Lieb-Thirring inequality. 

Lemma 3.3 (Lieb-Thirring inequality with positive Fermi level). Assume d > 1. 
Let < 7 < 1 be a self-adjoint operator on L 2 (R d ) such that |A + 1| 1/2 7|A + 1| 1/2 
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is trace-class. Then 7 is locally trace-class and its density satisfies 
(3.5) Tr |A + 1| 1/2 t|A + l| x / 2 > f 5T? c ( Pl {x)) dx 



where 



with po = \S d ~ 1 \ (2n) d /d, and where K(d) is a positive constant depending only on 
d. 



The proof of Lemma 13.31 follows ideas of Rumin [35 . Note that Lemma 13.31 is 
also valid in dimension d = 1. 

Proof. We follow a recent method of Rumin [35]. We introduce the spectral projec- 
tion P e := l(|A + 1| > e) in such a way that we have the layer cake representation 

/>oo 

|A + 1| = / P e de. 
Jo 

Let now < 7 < 1 be a smooth-enough finite rank operator. We have 

/>00 />00 n 

(3.6) Tr|A + l|7= / de Tr(P e7 P e ) = de p e (x) dx 

Jo Jo Jwi d 

where p e is the density of the finite-rank operator P e jP e . We now consider a 

bounded set 4cl' i and estimate 

. ,„ 2 

p e (x) dx = Tn{t A P e ^P e ) 



> 



W /2 



t A p e W /2 



(3.7) 



Note that, since II 7 II < 1, 

















62 / 



Tr(t A P e ^P e H A ) < HaP^L hi < \A\ /(e) 



where 



/(e) := (Shr)-* | { \ P 2 - 1 1 < e} | = ^ ((1 + e)^ - (1 - e)f ) 



d{2-K) d 

Taking A to be a ball of radius e — > centered at x, we obtain from (|3.7[) the 
pointwise estimate 

Pe{x) > (y/p^-Vm) 2 + . 
We may now insert this in (|3.6[) and obtain 

Tr|A + l| 7 > f dx I de(^fpjx)- ^/J{e)] 2 = [ R d (p(x)) dx 



with 
(3.8) 



Rd{ P ) 



de. 
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At zero we have 



Rd{ P ) 



(2n) d 2 



At infinity, one can compute that 



Rd(p) 



p^oo d + 4 



K sc (d) P 1+2 ' d . 



Hence there is a constant K(d) such that (|3.5[) holds. 

We have written the proof for a smooth enough finite-rank operator. The general 
case follows from an approximation argument based on Lemma l3.2l This completes 
the proof of Lemma 13.31 □ 

Since |A + 1| 1 / 2 Q ±± |A + 1| 1/2 is trace-class by assumption and < Q++ < 
II + < 1, — 1 < —II - < Q < 0, we immediately obtain from Lemma T3. 31 that 



(3.9) 



±Tr|A + 1| 1 / 2 Q ±± |A + II 1 / 2 > / Srr(p Q ±±{x))dx. 



It therefore remains to estimate the density arising from the off-diagonal terms Q + 
and Q h . 

Step 2. Estimate on Q ±T . It is enough to consider Q h = II - QII + , since Pq+- + 
Pq~+ — 23?pg-+. In order to estimate the density Pq-+ in the whole space K d , we 
argue by duality and write 



/ v PQ - + = Tr(i/n-Qn 4 



(3.10) 



Tr f — V - IA + 1I 1/4 0|A + ll 1/4 

llA + llV^ IA + 1IV4 + ^l A + 1 l 



n 



This calculation is valid if V is bounded and compactly supported, since Q is a 
smooth-enough finite-rank operator. Using Schwarz's inequality and that Q 2 < 
Q ++ — Q , we have 

"lA + ll^QIA + ll 1 / 4 _ < |A + 1| 1/2 Q 

©2 ©2 

1/2 



= (Tr|A + l|Q 2 ) 1/2 < (Tro(-A-l)Q) 



Returning to (|3.10j) . we obtain 

n+ 



v Pq - 

We now compute 



n+ _ y _ n- 



< 



v- 



ir 



IA + 1IV4 IA + 11 1 / 



IA + 1IV4 IA + 1IV4 



(2tt) 



-d 



(Tr (-A - l)Q) 



\V{p-q)\< 



1/2 



(3.11) 
where 
(3.12) 



(2tt) 



-d 



|Jg (i-H 2 ) 1/2 (W 

\V(k)\ 2 $ d (\k\)dk 
dp 

(l~b| 2 ) 1/2 (b-fc| 2 -l) 1/2 ' 



y^dpdq 



\p-k 



We will use the following fundamental result. 
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Lemma 3.4. For d > 2, the function $>d is bounded on R d . The function $1 is 
not bounded in a neighborhood of k = 2. 

For clarity the proof of Lemma 13.41 is postponed until the end of the proof of 
Theorem 12.11 We deduce that 

JaTW v \KTW ~ {2nrd/2 l<M ™ +) IMIl2( " 

which leads to the estimate 



e- 2 



(3.13) 



f \p Q - + | 2 < (2ir)- d |$ d || i oo (K+) Tr (-A - \)Q. 

JM d 



We can extend STf c for p < —po linearly as follows 

srr(p) = (po + P f+ ld (po) 1+ ~ d ^(PoFp- 

The function is now convex on the whole line R. Note that for d > 2, we have 
\p\ 2 > c5Ti c (p) for all p, hence we have also shown that 

(3.14) cf STT{p Q -+)<Tr (-A-l)Q 

Jm d 

for a small enough constant c > 0. 

Remark 3.2. Modifying the previous proof by using LI+| A + 1|~ Q V II~|A + 1| _Q 
with an appropriate power a and & p norms, one can show that 

(3.15) / \ PQ - + \ p +f \ PQ+ -\ P <C(d,p) Tr (-A-1)Q 

JR d JR d 

holds for all 2 < p < oo and all d > 2. 

Conclusion. Putting (|3.9p and (|3.14|) together, we deduce by convexity of STf c that 
Tr (-A - 1)Q > c / 57?° (p Q++ ) + <57? c (p Q -) + 57? c {p Q -+) + S7? c {p Q+ -) 

JR d 

> K{d) I srr ( PQ ) 

JR d 

for a small enough constant K(d) > 0. This completes the proof of Theorcm l2.ll □ 

Remark 3.3. Our method yields explicit values for the constant K(d) appearing 
in the statement of Theorem 12.11 see [7] . 



It remains to prove Lemma 13.41 

Proof of Lemma \3.4\ To study $>d(k) for d > 2, we make the decomposition p — 
(s,px) with s = p ■ k and find 

ds dp± 



(3.16) = |,S d - 2 | f 

J i . 



(s-kf+r^l (l- S 2-r2) 1/Z (( S -fc)2 +r 2_ 1 y 
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For k > 2 and |s| < 1, it is clear that (s — k) 2 + r 2 > 1. The integration domain is 
therefore independent of k when k > 2. It is then easy to verify that <£>d is decreasing 
and continuous on (2,oo). Hence we only have to prove that it is bounded in a 
neighborhood of [0, 2). Next we note that 



(8.171 <!>,,( A- 1 - \S d - 2 \ ( ) ,lr 



o VT^r* VVT 



r 



where we recall that 

dv r min(l-l+x) dy 



(v-xf>i VT^^iv-x) 2 - 1 J -i v 7 ! - ~rf^/{v-x) 2 - l' 



It is an exercise to verify that <&\ is a continuous function on R + \ {2} (in particular 
it has a finite limit at x = 0), and that 

*i(x) ~ —loglic-21, $i(x) ~ -. 

:r— ^2 z a:^oo X 

Using that, for instance, 

- , , c 

*i(z) < r > 

|x — 2| 4 



and letting u = y/i — r 2 , we obtain 

<E> (k)<C f 1 dU U <c(C dU Y ( f 1 dU 

" Jo VT^ 2 \k-2u\i ~ \Jo (l-u)i) \Jo \k-2u\i 

by Holder's inequality. The right side is bounded with respect to k, hence $ f j is 
uniformly bounded for d > 2. By a similar proof one can verify that &d is also a 
continuous function on R + . This completes the proof of Lemma 13.41 □ 



Remark 3.4. It is possible to calculate the exact maximum value of which 
might be interesting for physical applications [7 . Starting from (|3.16[) and letting 
t = r 2 , we obtain 

s , n f ds dt 

®3{k) = 7T / 2+t<1 — — — — 

J(s-kf+t>i {l-s 2 -t) 1/2 {{s-k) 2 +t-l) 1/2 

r>0 

min(l,fc/2) />1 — s 2 ^ 

ds 



mm(l,-l+fc) Jl-( s -ky (l - S 2 - t) 1/2 ((s - k) 2 + t - l) 
min(l,-l+fe) f 1 ^^ dt 

ds 



1 



(1 - S 2 - t) 1 ' 2 (( S - k) 2 + t - 1)V2 



In order to compute these integrals we use the fact that 



„ „n f dt a — t 

3.18 / ; = -2arcsin A / - 

J y/(a-t)(t-b) V a-b 

whenever a > b. We find for < k < 2, with s = — 1 + ku, 



$ 3 (fc) ^7r 2 + 2vrfc ^ ; 



/ 2 ku 7T 

( :!..!.<)) $ 3 (fc) =7r 2 + 27rfc ( / arcsin J u ^ - - du - - 
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We have 

2 - ku \ 2(u- 1) 



dk \2 + k[l- 2u) J (2 + k(l - 2u)) 2 

hence the function 



< 



t(k) := / arcsin \ u —. -du 

Jo V 2 + fc(l-2w) 4 

appearing in the parenthesis in p,19|) is decreasing with respect to k, by mono- 
tonicity of t arcsin(v^). Its value at k = is 



f 1 TT 
/(0) = / arcsin -\fudu — — = 0. 
Jo 4 



Therefore, /(fc) < for < k < 2. Now $ 3 (fc) = 2n(f(k) +kf(k)) < 0, hence $ 3 
is decreasing on [0,2]. Since we know already that $3 also decreases on [2, 00), we 
conclude that 

max $ 3 = $3(0) = 7r 2 . 

R+ 

Similarly as in (13.171) . we can express $^ in terms of $3 for d > 4 by assuming, 
for instance, k = k(l, 0, 0) and writing p — (q,p±) with 96I 3 and p± £ M. d ~ 3 . 
We obtain the recursion relation 

(3.20) $ d (fc) = IS"'" 4 ! / y/l -r 2 $ 3 ( k \ r d ~ 4 dr, for d > 4. 

Jo Vvl-r 2 / 

As we have shown above that $ 3 is strictly decreasing, this proves that $d is also 
strictly decreasing, hence that 

max$ d = $ rf (0) = 7r 2 |5 d " 4 | / \/l - r 2 r d ~ 4 dr, for d > 4. 
K+ JO 

3.3. Proof of Theorem 12. 2t kinetic Lieb-Thirring inequality for d = 1. In 

the one-dimensional case d = 1 , the same proof as that of Theorem 12.11 leads to a 
bound of the form 

(3.21) Tr (-A - 1)Q > K(l) [ 6T C (p Q++ + p Q — ) 

JR 



+k{1) Jr *m — 



|2 



Using the known behavior of $1 at \k\ = 2 and when |fc| — > 00, one can state this 
bound as in ([2~T4]) . □ 



4. Potential inequalities: Proof of Theorems 12.31 |2~H and |2~51 

For the standard Lieb-Thirring inequalities [26, 27] (the case where fi = 0), there 
is a duality between the kinetic energy and the potential versions of the inequality, 
and this duality is based on a variational principle for sums of eigenvalues. A similar 
variational principle is also valid inside the continuous spectrum and can be used 
to deduce Theorems |2~B1 and |2~41 from Theorems |2~T1 and |2~2"1 
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Theorem 4.1 (Variational principle). Let p > and V be a real-valued function. 
Assume that V G L 2 (R d ) n L 1+d / 2 (R d ) when d>2, and that V € L 3 / 2 (R) + L 2 (R) 
with 

whend = 1. Then bothQy \ — A—p] 1 / 2 andQy | — A — p+V\ x ^ 2 are Hilbert- Schmidt 
operators and 

(4.2) Tr v (-A-n + V)Q v = inf ( Tr (-A - /i)Q + / y(a;)pQ(x)dx 

where K, was defined in Q3.3p . TTie infimum in (|4.2[) is attained for Q = Qv ■ 

To motivate this theorem, we explain its analogue for self-adjoint finite-dimensional 
matrices A and B. The starting point is the well-known formula for the sum of 
eigenvalues [2TJ Thm. 12.1] 

-Tr(A + S)_= inf Tr(A + B)-y. 

0<7<1 

Introducing the spectral projection II = 1(A < 0) onto the negative spectral 
subspace of A and changing variables, 7 = Q + II - , we obtain 

- Tr(A + B)- = inf Tr(A + B)Q + Tr(A + B)IT , 
-n-<Q<i-n- 

that is, with the notation IT^ = 1(A + B < 0), 

Tr(A + B)(np -IT) = inf Tr(A + B)Q. 

-n-<Q<i-n- 

The right-side is obviously the analogue of the corresponding term in (|4.2|) . with 
A = —A — p and S = V. The left-side is negative, which can be seen by taking 
Q = on the right, or by noticing that 

Tr(A + B)(n B - IT) = - Tr \A + B\(Ug - II-) 2 . 

This, clearly, is the analogue of Try(— A + V — p)Qv, see Definition 12.21 



4.1. Proof of the Lieb-Thirring inequalities in a potential V. Here we ex- 
plain how to prove the Lieb-Thirring inequalities (|2.24[) and (|2.28[) . assuming The- 
orem [47T] As in the proof of Theorem 12. 11 we assume p = 1, the general case being 
obtained by a simple scaling argument. By Theorems 14.11 and 12 . 1 1 we have for d > 2 

> Try(-A + V — l)Q v 
> inf (K{d) [ 6T SC { P )+ f Vp) 

= -L(d) {{V{x) - 1)^ - 1 + ^ V(x)) dx. 



The second equality follows from a simple optimization argument. When d = 1, 
we argue similarly. We decompose p = pq++ + Pq— and p' = Pq- + + Pq-+ and 
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use (13.21[) to obtain 

> Tr y (-A + V - 1)Q V 



> inf (k(1)[ ST sc ( P )+ I Vp)+M[K'(l)[ l p^dk+ [ V p' 
p>-po \ Jr Jr ) p ' \ Jr Fi{\k\) Jr j 

= -L(l) / ((V(x) -l)f-l + | V(x)) dx - L'(l) [ ^(|fc|)|\>(fc)| 2 dk 



with 



and 



F i{\ k \) = ,,, log 



\k\ °\\2-\k 



3 \3K(l) J ' w 4JT'(1))' 

This concludes the proof of (pEMj) and (ET2"51) . □ 

4.2. Proof of Theorem I4.lt the variational principle. As before we assume 
/i = 1. Let us denote by I(V) the infimum appearing in (|4.2[) : 



I(V) := ^( Tr (-A-1)Q+ / 7(a;)/jQ(x)dx 



Note that by Lemma I3T21 we can restrict the infimum to finite-rank states Q e JC. 
We split the proof of the theorem into two parts. First we show that 

(4.3) > Tr y (-A - 1 + V)Q V = - Q v \ - A - 1 + V\ 1/2 2 > I(V) 

This will show that Q v \ - A - 1 + V\ 1/2 is Hilbert-Schmidt. We will also find 
that Q v \ - A - ll 1 / 2 is Hilbert-Schmidt, To prove (|4.3p , we approximate Qy 
by a well-chosen sequence Qy of smooth operators in X satisfying the constraint 

-n- < Q v < n+. 

In a second step we prove the converse inequality 

(4.4) Tr v (-A-l + V)Qv <I(V), 

using the information that Qy \ — A — 1 + V] 1 ^ 2 € 62 and the density of finite-rank 
operators in IC, as stated in Lemma I3T21 

Step 1. Proof of the lower bound (14. 3|) . We introduce the following function 

h(x) := -\x\ t(\x\ < 1) + (2 - \x\) 1(1 < |x| < 2) 

and replace —A by K e := —A + eh(-iV) for a small e > 0. The gain is that 
—A + e h(—iSJ) now has a gap (1 — e, 1 + e) in its spectrum. Note also that we have 
IT~ = t(K E < 1) for all e > 0, hence the free Fermi sea is not changed. 
Let us introduce the corresponding regularized operator 

Q e v := l(- A + e/i(-iV) + y < l) -IT. 

Note that — A + eh(— iV) + V — > — A + V in the norm resolvent sense. When d > 3, 
our assumption that V € L 2 (R d ) n L 1+d / 2 (R d ) implies V € L( d+1 )/ 2 (M d ) hence it 
follows from a result of Koch and Tataru [15] that — A+V has no positive eigenvalue. 
This in turn implies that Q v -> Qy strongly by e.g., [32, Thm. VIII. 24]. 

In dimensions d — 1 and d = 2, it was shown by von Neumann- Wigner 39] and 
Ionescu-Jerison [JJ] that there exist potentials V satisfying our assumptions for 
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which kcr(— A + V — 1) ^ {0}. For this reason, when d = 1, 2, we assume first that 
V has a compact support and is bounded (then — A + V has no positive eigenvalue 
and Qy — > Qy strongly), and only remove this assumption at the very end of the 
proof. 

What we have gained is that the operator Q v is now Hilbert- Schmidt, whereas 
Qy is not even compact in general (Remark 12. 4[) . 

Lemma 4.1. Under our assumptions on V, 

(4.5) Qye© 2 , (-A)Q^e6 2 and (K E + V) Q £ v e 6 2 
for all e > 0. 

For clarity we postpone the proof of Lemma 14.11 

The facts that Q v € 6 2 and (-A)Qy £ 6 2 imply in particular (Qy) ±:t € ©i, 
|1 + A\ 1/2 Q e v E 6 2 and |1 + A| 1 /2(Qe r )±±| 1 + A |i/2 £ 6i _ In particular Q £ v e X, 
the Banach space introduced before in (|3.1j) . By Theorem 12. II in dimensions d > 2, 
we deduce that pq^ is a well-defined function such that STi c (pq v ) € L 1 (R d ) and 
that 

Tr (-A - l)Q v = Q v \ - A - lf^ 2 > f ST?^). 

In dimension d = 1 we at least know that pq<s € Lj L oc (M) by Lemma 13. II 
Lemma 4.2. We have the following equality: 

(4.6) Tr (-A- l)Q v +e Ty\h(^)\((Q v ) ++ - (QV)-') + [ VpQ v 

= - Tr | K e + V - 1 1 1/2 {Q v f \K e + V - 1 1 1/2 < 0. 

Proof. It is possible to approximate Qy by a sequence {R n } of smooth finite rank 
operators such that —II - < R n < U + , (—A + i)R n —> (—A + i)Qy strongly in 6 2 
and (-A + i){R n ) ±:t (-A +i) -> (-A + i)(Q y ) ±=t (-A + i) strongly in 6i. See, 
e -g-> LSI Prop. 2 & App. B]. We then write 

Tr(— A - 1 + eft(-iV) + V)i?„ = Tr(— A - l)R n + e Tr h(-iS7)R n + / F/j fln 

= Tr |tf e - 1 + V\ (n+ £ i?,„n+ B - Uy e R n U V£ ) 
where IIy £ = t(K £ + V < 1) and pass to the limit n — > oo. □ 
Since Q v € <-f, we deduce from (|4.6p that 

(4.7) - Tr \K e + V - \\ 1/2 {Q v ) 2 \K e + V- 1| 1/2 > /(V) 

for all e > 0. In particular Q v \K S + V — l] 1 / 2 is uniformly bounded in the Hilbcrt- 
Schmidt class @ 2 . Note that the weak limit of Q V \K £ + V — l] 1 ^ 2 in <3 2 can only be 
Q v \-A+V-1\^ 2 , since Q e v ->■ Qy strongly and 1^+^-1 1 1 / 2 ?/ ™> l-A+F-ll 1 / 2 ?/ 
for every y g i? 2 (IR d ). This latter statement can be seen by writing^ 

lK , + „ _ ir „ . iy;-'' 1/a ( -A + v + » +e | y + Vi |1 ^ Hv ) » 

~ A + 1 ' 1 1 ' 2 (-A + V + i)y = | - A + V - l\^ 2 y 



e^O -A + V + i 



In dimension d = 1, the domain of — A + V (hence of K E + V) contains H 2 {K) by choice of 
the Friedrichs extension via the KLMN Theorem. 
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and using that f(x) = \x— l| 1 / 2 (a; + i) 1 is a continuous function tending to zero at 
infinity, thus \f{K e + V - 1) - /(-A - 1 + V)\\ -> by [H2 Thm. VIII.20]. Hence 
we have 



^-H 1/2 Qf/ 



A + y-l| 1/2 Q y weakly in 6 2 



and, passing to the weak limit in (|4.7[) . we obtain the claimed inequality (|4.3 
From (|4.6[) . we also have the following bound 



(4.8) 



IA + 1 



1/2 



J V PQ% , = Tr (-A - l)Q v + J Vp Q s v < 



for all e > 0. We deduce for instance that 



(4.9) 



Qv|A + l 



1/2 



< 



Tr (-A - l)Q v + 2 y Vp Q ^j < -I(2V). 



This uniform bound proves that Qv\A + ll 1 / 2 G ©2 and that 

Q^|A + 1| 1/2 Qy|A + II 1 / 2 weakly in & 2 . 

In dimensions d — 1,2, we have only written the proof for V a bounded function 
of compact support. If V is an arbitrary function satisfying our assumptions (I2.16[) 
and (|2Tf|) . we apply the result to V R (x) := V(x)l(|x| < R)t(\V(x)\ < R) and, 
from (|4.3p and (|4.9|) . we obtain uniform estimates of the form 



and 



Tr I - A + V R - l\ 1/2 (Qv R ) 2 \ - A + V R - l^' 2 < -I(V R ) 



Tr |A + 1| 1/2 (Q++ - Qy~)\A + l^ 2 < -I(2V R ). 

Extracting subsequences we now have at best that Qv R ~ ^ weakly as R — > 00, 

where < S < 1(— A+V = 1). Passing to weak limits as before, we therefore obtain 
that 

Tr I - A + V - l\ 1/2 (Qv) 2 \ -A + V- 1| 1/2 < -I(V) 
as was claimed, and that 

Tr |A + 1| 1/2 (Q++ - Qy~)\A + \\ x ' 2 < Tr |A + 1| 1/2 (Q++ - Q v ~ + S)\A + \\ x ' 2 

< -I(2V). 



It remains to provide the 

Proof of Lemma \4-l\ Our claim (|4. 5|) follows from Cauchy's formula and the resol- 
vent expansion: 




Under our assumptions the function V is -fQ-compact, hence K e + V has the gap 
(1 — e, 1 + e) in its essential spectrum and it is bounded from below. In (|4.10[) . 
we choose for C a smooth curve enclosing the spectra of K e and K e + V below 1, 
without intersecting them. We will explain below how to choose J. 
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In order to show that (1 — A)Q V is a Hilbert-Schmidt operator for all e > 0, we 
estimate each term in (|4.10p . Our bounds will depend on e. We start by noticing 
that there is a uniform bound of the form 



(4.11) 



Vz£C, 



1- A 




K E -z 


+ 



K £ + V 



< C. 



The constant C diverges when e — > but we do not emphasize this in our notation. 
To estimate the last term of (I4.10[) . we use (for d > 2) that 



(4.12) 



V- 



K F -z 



<C\ 



©l + d/2 



\\L 1 + d / 2 ( 



by the Kato-Seiler-Simon inequality [36 , Thm 4.1], 



(4.13) Vp > 2, ||/(-iV)ff(a:)|| 6 , < — tIML^H/IIl^)- 

(2ir) p 

The constant C in (|4.12j) also depends on s. Choosing J > 1/2 + d/4 in (|4.10|) . we 
obtain by Holder's inequality and (|4. 



(1-A) 



V- 



J+i 



{Ke - Z, 



1 



K E + V 



■ dz 



<c\\v\\ J L t^ /2( 



We now treat the term corresponding to k = 1 in the first sum of (|4.10l) and start 
by noticing that 



IT 



-V- 



IT 



z K F 



dz = 



^^ 



^^ 



—V- 

z K F 



■ dz = 0. 



For the other terms, we simply write, for instance, 



(1-A) 



IT 



-V- 



^^ 



K F - z K F -z 



dz 



<c\\urv\\ 



<c\v\\ m 



since II - = 1(|V| < 1) has a compact support in Fourier space. 

The argument is the same for the other terms of the first sum in (|4.10|) : We 
write 



IT 



IT 



K £ — z K £ — z K £ — z 

and note first that the term which has only Il + vanishes after integrating over zeC, 
by the residue formula (the same holds for the term which has only IT~ ) . The other 
terms contain at least one II - and can be estimated similarly as before. 

We deduce, as was claimed in (|4.5j) . that (1 — A)Qy £ ©2 for every e > 0. Since 
{Qv? = (Qv) ++ -(Qvy~> this implies that (1- A)(Qf,) ±=t (l- A) £ Si. Finally, 
(K £ + V)(l - A)" 1 being bounded, we have that (K £ + V)Q V £ 6 2 - 

We have written the proof for d > 2. The case d = 1 is similar and left to the 
reader (see also the proof of Theorem 12.51 below). □ 
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Step 2. Proof of the upper bound (|4.4j) . To finish the proof, it remains to show the 
inequality flO}, that is Try(-A - 1 + V)Q V < I(V). 

We pick a smooth finite rank operator Q such that —II - < Q < H + and Q{— A) 
is bounded, and note that 

Tr (-A-1)Q+ f V PQ 

= Tr(— A - 1 + V)Q 

= Tr | - A - 1 + y| 1 / 2 (n+Qn+ - U V QU V )\ - A - 1 + V\ 1/2 . 
We now use that 

| - A - 1 + V\ 1/2 U V Q V U V \ - A - 1 + V\ 1/2 e 6i 
as we have shown in Step 1. Writing Q = (Q — Qv) + Qv we obtain 

Tr (-A-1)Q+ f Vp Q = Tr v (-A-l + V)(Q-Qv) + Tr v (-A-l + V)Qv 

>Tr v (-A-l + V)Qv 
In the second line we have used that 

Trv(-A-1 + V)(Q-Q v ) >0 

since — Ily < Q — Qy < Hy. By the density of finite rank operators in X (see 
Lemma I5T^|) , we deduce that 



Tr v (-A + V-l)Qv<I(V), 
which finishes the proof of Theorem 14.11 
Remark 4.1. Our proof also yields the limit 



□ 



(4.14) 



lim 

£^0 



K e + V - 1 



1/2 



Q e v 



V -1 



1/2, 



Indeed, from (|4.7[) . we know that 
lim sup \K £ + V- l| 1/2 Qy 

e-i-0 



<-i(v)= I 



V-l 



1/2/ 



where the last equality follows from Theorems l2.3l and l2.4l Since we also have proved 



that \K £ + V - 
follows. 



1|V2< 



A + V— l\ 1/2 Qv weakly in 6 2 , the statement (|4TT]) 



4.3. Proof of Theorem 12. 5t second-order perturbation theory. In this sec- 
tion we sketch the proof of Theorem 12.51 We detail first the one-dimensional case 
d = 1 and mention the necessary modifications in higher dimensions afterwards. 

We could embark upon expanding Q t v in powers of t by directly using the 
resolvent formula. Since we want to avoid a tedious justification of this expansion, 
we instead work with the approximate state 



itv 



tV < 1) -LT 



which we have already introduced in the proof of Theorems 12.31 and 12.41 We will 
prove bounds in t which are uniform in e, and pass to the limit e — > in the end, 
using (|4.14p . The same method of proof can be used to justify an expansion of Qtv 
to any order in t. 
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We come back to the resolvent expansion (|4.10p for Q\ v which we have already 
mentioned in the proof of Lemma 14. 1 1 above. In dimension d = 1, we write 

(4.15) 



Qlv =*Q! + t2 Q E 2 + t 3 Ql(t) 



where 



1 



2l7T J c K, 



-V- 



K, 



■ dz, 



and 



K K 



Ql = -77— 



-V 



1 

2l7T J C 
1 



K K 



-V 



■ dz 



K e +tV 



■ dz. 



In the above formulas, we choose for C a curve in the complex plane enclosing the 
interval [-R, 1] C E, where -R < inf a(K E + tV) for all < e < 1 and all \t\ < 1. 
To simplify certain estimates below, we also assume that \^sz\ < 1/2 for all z E C 
(in such a way that log > 0). For convenience we will make the assumption 

that 1 ^ <t(K s + tV) for all t small enough. If 1 is an eigenvalue of K £ + tV, one 
has to let the curve C depend on e, and modify it a bit in a neighborhood of z = 1. 
It can then be verified that our estimates below still hold true. These details are 
left to the reader for brevity. 

Note that Q\ is purely off-diagonal, i.e. (Q\ ) ±±: = 0. Using that \K e — z\ > for 
all z 6 C by definition of K s , one can prove (similarly as in the proof of Lcmma [4.ip . 
that Q\, Q\ and Q%{t) are trace-class, and that 

(4.16) Tr tv {K E - 1 + tV)Ql v = t 2 { Tr(QfU) + TrQ e 2 (K s - 1)} 

+ t 3 { Tr(Q e 2 V) + Tr Q%{t){K s -1+tV)}. 

Each of the terms of the right side makes sense and can be bounded uniformly in t 
and s, as we now explain. First, we have 



\\q £ 2v\\ 6i < (^r 1 f c 



i 



K F - z 



-V 



\dz\, 



6i 



and, similarly, 



\\Qt{t){K s -l+tV)\\ ei <C j 



K F -z 



-V 



\dz\, 



Si 



since (K e + tV — l)(K e + tV — z) 1 is uniformly bounded for z £ C, by choice of 
the curve C in the complex plane. We now use that 



A- II 



< C 



K e - Z 

for a constant C independent of z € C and e, to deduce that 



K F — z 



V 



<C\%z\- 1/2 



1 

- A-Z|V2 



v- 



1 



11/2 



- A-zP/2 



V 
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We have the bound 

1 1 
:V 



A-z] 1 / 2 \-A-z\V 2 



(4.17) 

and, in a similar fashion, 



= (27T)- 1 ||F|| il(E) | r 
<C\\V\\ Ll(m hg\Zz\-\ 



dp 



A-z] 1 /' 



■V 



<C\\V\\ L2m log\Zz 



Using these two bounds we deduce that 



Kr-Z 



V 



<C\V\\ im \V\ mm ^(-^(loglSzr 1 ) 



Si 



Integrating over zgC, this eventually shows that 

\\QW\\ 6l + \Ql{t){Ke - 1 + tV)\\ 6i < C \V\\ Hm) \\V\\ L2{m , 

hence that 

(4.18) | Tt w {K e - 1 + tV)Ql v - t 2 { Tr(Q\V) + TiQ e 2 {K e - 1)} 

+ 3 || T/ ||2 



<ct A \v\\< \\v\\ LH 



with a constant C independent of e and t. 
Using the residue formula we find 



Tt(QIV) = -2(2tt)- 1 



\p\'<1 \q\ 2 -\p\ 2 +e(h(q)-h(p)) 

I<?| 2 >1 



dpdq 



and 



TrQ E 2 (K c -l) = (2n) 



\V( P ~q)\ 2 



\P\1<1 \ q \2-\ p \2 + s(h(q)-h(p)) 



dpdq. 



The result in the case d = 1 now follows from taking first the limit e — > in (|4.18l) . 
using (|4.14|) . and then t 0. 

When d > 2, the proof is similar but a bit more tedious. We start again with the 
resolvent expansion (|4.10l) , to an order J such that the last term becomes trace-class 



30 RUPERT L. FRANK, MATHIEU LEWIN, ELLIOTT H. LIEB, AND ROBERT SEIRINGER 



(when multiplied by K e +tV— 1). This means we write 



Tr ty Q E tv (K e +tV-l)=- t 2 (2ir)~ d 



\V(p-q)\ 2 



\P\*<1 \q\2-\p\ 2 + s(h(q)-h(p)) 

\q\ 2 >l 



dpdq 



J— 



-v 



-V 



K e -l 

K F — z 



(4.19) 



3=3 

+ (-t) J+1 Tr j> dz 



1 



-V 



K E -z 

We fix a J > 1 + d/2 and deduce, similarly as before, that 

©i 

i 



J+i 



K e + tV — 1 



if £ -z ) 
< C|3z|" 1/2 



1 



j 

Sl + d/2 



-A-*|Va |_A-«|Va 

x (\\Vl L * m log\%z\- l + \\V\\ L 



A - zP/2 



i<*) 



with a constant C that is independent of e. For the other terms in (|4. 19[) . we have 
to work a bit more. As an illustration, we only consider the term 



Ti n <t> dz 



1 



K P -z 



-V 



the other terms are treated by the same argument. We decompose 

i n- n+ 



K F - z K e -z K s -z 

and expand {{K e — z)~ 1 V) 3 accordingly. The terms which have only II + or only 
II - vanish after the integration over the curve C, by the residue formula. For the 
other terms, n _ VTI + (or its adjoint) must appear at least twice in the trace to be 
estimated. For instance, we look at the term 
(4.20) 

n+ n+ n- / n+ n+ n 

Tr d> dz V V V = Tr & dz— V— V— VIL+. 



K F - z K F -z K F 



K F -z K P — z K F 



By cyclicity of the trace, this term can be estimated by 
(4.21) 



I62DH< jf \ d A 



n+ 



:V 



^^ 



A-zP/2 l-A-zP/2 



n+ 

- A-«|V2 



V- 



rr 



11/2 
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Decomposing \-A-z\^ 2 = | - A- z| 1 / 2 l(|A + 1| > 1) + 1 - A- z| 1/2 l(| A + 1| < 1) 
and using that V £ L 1 (M d ) n L 1+d ^(R d ), we find 



ir 



A-zl 1 / 2 
< 



I'- 



ll 



z |l/2 



1 



< 



- A-*|Va 
IdA + ll < 1 



11/2 



e 2 



1(|A + 1|>1) 



11/2 



< C (lVl L1(Md) log \Qz\- 1 + \\V\\ L1+d/2(Md) j . 
For the second term in the right side of (|4.21|) . we use that 

n+ „ n 



z |i/a 



V- 



z\V 2 



< 



n+i(|A + i| < i] 



11/2 



V 7 



rr 



A 



11/2 



n+i(|A + i| > l) 



A 



11/2 



V- 



n 



11/2 



The first term on the right side is estimated as before. For the second one, we use 
that 



n+l(|A + l| > 1) 



- A-zP/2 



V- 



n- 



- A-zP/2 



< C\%z 



n+ 



i] 1 /- 



-v- 



n- 



1|V4 



This term is now exactly the one which we have calculated before in (|3.11|) and 
it is finite under our assumptions on V. Summarizing, we have proved that the 
term (|4.20[) is bounded uniformly in e. 

The same argument can be applied to all the terms in f|4. 19[) , showing that they 
are bounded uniformly in e. This concludes our sketch of the proof of Theorem l2.5l 

□ 



5. Thermodynamic limit and positive temperature 

5.1. Lieb-Thirring inequalities in a box. In this section, we extend our in- 
equalities (|2.9[) and (I2.24p to the case of a system living in a box of size L, with 
constants independent of L. For simplicity we restrict ourselves to periodic bound- 
ary conditions and dimensions d > 2. 

We denote by —Al the Laplacian on Cl = [— L/2, L/2) d , with periodic boundary 
conditions, and, for any chosen \l > 0, we introduce IT£ „ := 1(— Aj, < /i). Note 
that since the spectrum of — A^ is discrete in has finite rank for every finite 

L > and fi > 0. The following is a generalization of the density inequality (12.91) . 

Theorem 5.1 (Lieb-Thirring inequality in a box, density version, d > 2). We 
assume that d > 2, /j, > and L > 0. Let Q be a self-adjoint operator of finite rank 
such that — LT£ < Q < 1 — . Then there exists positive constants K(d) and C 
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( depending only on d > 2) such that 
(5.1) Tt L 2 {Cl} (-A l -^)Q 



> K(d) 



where we recall that 




Ipq^I-CL- 1 ^) ) dx for p> 1/L 2 



<5T SC ((|pq(x)| - CL- d ) + ) dx for p < 1/L 2 , 



5V C ( P ) := (po + P ) 1+1 d (po) 1+1 * ^(Pof d P 

with po = p^ q{2Ti)- d \S d - 1 \/d. 

The function appearing in the integrand of (I5.1[) vanishes for p < CL~ 1 p~^~ 
(in the case p > 1/L 2 ) or for p < CL~ d (in the case p < 1/L 2 ), and it converges 
to STh(\pq\) in the limit L — s- oo. Note the absolute value which we have used to 
simplify our statement. Of course, ST^ c (\pq\) is comparable to 5T^ c (pq). 

Using Theorem 15.11 we can now deduce the (dual) potential version in the box. 
Again, note that for V <E L 1+i1 / 2 (Cl), the spectrum of — A^ + V is discrete and 
bounded from below, hence there is only a finite number of eigenvalues below each 
chosen Fermi level p. 

Theorem 5.2 (Lieb-Thirring inequality in a box, potential version, d > 2). Assume 
that p > 0, d > 2 and L > 0. Let V be a real-valued function in L 1+c1 / 2 (Cl)- Then 
we have 



(5.2) > - Tr(-A L + V - p)- + Tr(-A L - M )_ - p [ 

Jc 



V 

d-l 



-LUI< j J (V(x)~pt + ^ -p 1 + -2+^piv( X ) + ^\Vir) 1 d.r 



2 ^ v ' L 



when p > 1/L 2 , and 



(5.3) > - Tr(-A L + V - p)- + Tr(-A L - - p [ V 

Jc L 

> -L(d) (v(xf^ + j- d \V{x)\^ dx 

when p < 1/L 2 . The constant L{d) only depends on d. 

Since all operators are finite-rank, the proof simply reduces to computing the 
Legendre transform of p H> (l + {\p\ — e) + ) — 1 — a(\p\ — s) + . We skip the details 
and only provide the proof of Theorem 15.11 



Proof of Theorem \5.1\ The proof follows the same two steps as that of Theorem l2.11 
but it is slightly more tedious. 
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Step 1. Estimate on Q . We start by estimating the diagonal densities Pq±±- 
Following the strategy of the proof of Lemma \3. 31 we get, with 7 = Q , 

Tr i 2 (c L ) I - A L - /1I7 > / Rd,n,L(p-,( x )) dx 
Jc L 

where 

(5-4) Rd^ L (p)=J o (yp~-\]fd,»,L{e)) de 

and 

r ( e \ - 

L d 



W(e) = -L # {p G (2nZ/L) d : \p 2 - /x| < e} 



1 

~ l7 1 ^ ' ' " \ 2vr V" ' ^ 7 I " " ' " I 2tt V" A 4 

The following gives an estimate on the function fd,n,L- 

Lemma 5.1 (Estimates on fd.^.h)- When /1 > l/L 2 , we have 

I d-i \ 
(5.5) /d,M,i(e) < C i^- + M l_1 e l(e < /x) + e d / 2 l(e > /i) 



whereas when < /1 < l/L 2 , we have 

r.l 

L d 



(5-6) / <JlM , £ (e)<C(-l + e''/ 2 



/or all e > 0. 

Note that the estimate (|5.6I) on fd.^.L in the case /1 < l/L 2 is a bit weaker than 
the one (f575l) for /i > l/L 2 . 



Proof of Lemma \5.1\ First, we recall the following well-known property 



(5-7) 



#Z d n B(R) - —r^R' 



< Cmax(l,R d ~ 1 ), 



which says that the number of points of the lattice 7L d inside a ball of radius R, 
behaves like the volume of the ball B(R) in the limit of large R, whereas it is just 
bounded for small R. The error term can even be replaced by o(R d ^ 1 ) but we do 
not need this here. Note that the bound ([ETfjl implies #Z d n B(R) < C(l + R d ). 
The proof of ()5.6() is now straightforward: Assuming [i < l/L 2 , we simply write 

We) < ^* Zd B ih^) ^T d ( 1 + Ld ^ 12 + e " /2 )) ^ I? + Ced/2 - 



In order to prove (|5.5[) we need another estimate. Let M > and < x < Xq 
for some fixed xq > 0. Using (|5.7I) we obtain 



#z d n b (m(i + x) 1 ' 2 ^ -#z, d nB (m(i - x) 1 / 2 ^ 

S " " ((l+x) 1 / 2 -(l-x) 1 f)+2Cm a x(l,M d - 1 )(l + x y 



d 

IS" 1-1 ! M d 



< C- j x + 2Cmax (1, M^ 1 ) (1 + x ) ! 

d ' 

(5.8) < C(M d x + max (l,M d_1 )) . 
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We have used that (1 + a;) 1 / 2 — (1 — x) 1 ^ 2 < Cx for all < x < xq, where C only 
depends on x$. We can use (|5.8[) to prove (|5.5p . assuming now p > 1/L 2 . For 
e < 3^/2, we use (JEHJl with M = L^/JI/2ir > 1/(2tt) and x = e/fjt < 3/2. We 
obtain 

W«) < ^ + ^V-j < C + L j . 

Finally, for e > 3p/2 we have 

l/2\ 



fd, 



< C ( M d/2 + e d/2 + ) < C | — + e 



d-1 

1 \ ^ A 1 " 



L d J ~ \ L 

where in the last estimate we have used both L -1 < p 1 ! 2 and e > 3/i/2. This 
finishes the proof of Lemma 15.11 □ 

Using the bounds (|5.5[) and (|5 . 6[) on fd,fj,,L, we can now deduce an estimate on 
Rd.^.L appearing in (|5.4p . To simplify our argument, we introduce 

f^K^e^rf + e^lfe^/i) for ^ > 1/L 2 , 
(5 - 9) ^^-{e^ for M < 1/L 2 . 

such that (|5.5p and (|5.6I) can be rewritten as 

/d, M ,L(e) < £d, M ,L + Cg d ,^,L(e) 

with 

L~ 1 p i i 1 for /! > 1/L 2 , 



We then have in all cases 



L~ d for /i < 1/L 2 . 



Rd^Ap) = (yp- \J fd.^d e ) s j de 



with 



Sd^.L(p) = P - \ 9d,^,L{e) de 



>C 



'p^p 2 t(p < p 2 ' d ) + p 1+2 / d l(p > p 2 ' d ) for /i > 1/L 2 , 

p 1+2 / d for p< 1/L 2 . 

To conclude, it suffices to note that 

{VP- yf£djl) 2 + > a g (p~0e difliL ) + 

for any 6 bounded away from and a$ small enough, and that Sd^^hiptep) > 
Pe S d ^x{p). 
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Step 2. Estimate on Q ±T . We again separate the cases p > 1/L 2 and p < 1/L 2 . 
We start with the case p > 1/L 2 and decompose Q H as 

Q+- = n+Qir = (n+ + n+)Q(n - + nr) = QJr + Qtr + Qto~ + Qtr 



where 

n„ = 1 (p < P 2 < P + vW X ) and n o~ = 1 (m ^ vW^ <P 2 < p) 
(we remove the index on IF L for simplicity). We have, with := £ _d / 2 e lfe ' :E ; 

fe,£e(27rZ/L) d 
/i<fc 2 <^+ v /7T/L 
^i-^/L^f 2 ^ 

The matrix (efe.Qe^) has a norm < 1, hence we deduce by Schwarz's inequality 

that 

(5.10) 



1 



J 1 



In the last bound we have used (|5.8|) and the assumption that p > 1/L 2 . For p^ - , 
we write, this time, 



! Tr(VQ 



10 



A L -/i| 1/2 Q 



< v/Tr(-A L - M )Q 



1/2 



We now have 



^l-Ai-HVa 



1 



inp-g)l 2 

g 2 >/j+v^t/t 



E 



< 



1 L 



E \v{k)\ 2 <c^ f \v\ 



|p-fc| 2 >At+VP/i 



where, in the last estimate we have again used that 



L~ d #{p-^/L< P 2 <p) <C l - 



From these bounds we deduce that 



(5-11) 



C L 



pt -\ 2 <Cfi—Tr(-A L -v)Q- 



The term p^i is treated similarly. We conclude this paragraph with an estimate on 
pl^ , which we derive by the same method as for (|3.13j) . in the proof of Theorem l2.ll 



(5.12) 



c L 



d) Tr L 2 (Ci) (-A L - l)Q, 
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with 

(5 ' 13) : =ET £ i 1^1/^1 , |2 ^ - 

P e(2irz/ L y (^-bl 2 ) (b-fc| 2 -/i) 

|p-/c| 2 >Ai+VM/L 

The function <&d.pi,L is a Riemann approximation of 2 )/ 2 <I> d (- / In order to 
prove that $d,/x,L is uniformly bounded on (27rZ/L) d by C/i( d ~ 2 ^ 2 , independently 
of L, we compare it with its limit. For every p in the sum above, we introduce the 
ball Bp of radius r]/L, centered at p. We will fix the value of ij later, but as a first 
constraint we impose that 



(5.14, V „-f + a„ d ^ + f-f^ + f. 

for all /i > 1/L 2 and L > 1. It is easy to verify that the previous condition is 
satisfied when, for instance, rj < 1/8. The constraints f|5 . 14[) imply that 

(5.15) V P >eB p , \ p f<„-^-, \ p >- k f 

Next we compute the gradient 
1 



%-bl 2 ) 1/2 (b-fc| 2 -^) 1/2 

p p — k 



— If I 2 b-fc| 2 -M/ (^-bl 2 ) 1/2 (b-fc| 2 -^) 1/2 ' 

For satisfying (|5.15[) , we have 



/.-|r>T~ YE v 2i %^ 



and 



|p'-fc| 2 -^- ^ V 2L ^m 

We therefore deduce by Taylor's formula, that for every p' G B p 



(m- bl 2 ) 1/2 (b-fc| 2 -M) 1/2 (M-b'l 2 ) 1/2 (b'-fc| 2 -M) 1/2 

< (2 + V6) (2r?) sup 1 



Choosing 77 small enough, we can therefore make sure that 



1 1 



SU P — T7o : -T-T < 2 



?eB p (^-IgP) 1 / 2 (l^fcp-^ 1 / 2 (/,-b| 2 ) 1/2 (Ip-fcp-^) 1 / 2 

and then that 

1 



< 4 inf 



(/* - bl 2 ) 1/2 (b - fcl 2 - m) 1/2 (/* - M 2 ) 1/2 (k - ^l 2 - m) 1/2 ' 
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Using that the balls B p are disjoint for rj small enough, we finally obtain 



4 



E 



i 



bl 2 <^-vW-£* 

\p-k\ 2 >H+y/ji/L 

^ \ S d-x\rfi ^®d{\k\/yfji) < Cp^, 

since $>d is bounded by Lemma l3.4l Summarizing all our estimates, we have proved 
that 



m d P' 



pt -\ 2 <CTr L , (CL) (-A L -ii)Q. 



Using now both that \x — e\ > (\x\ — e) + and p$ < Cp( d x ^ 2 /L, we deduce that 



/' 



c L 



CvP-W/L) <CTt l2(Cl) (-A l -h)Q 



with a constant C that does not depend on L, for p > 1/L 2 . This completes the 
proof of (|5.1j) when p > 1/L 2 . 

The case p < 1/L 2 is similar, except that we only decompose 



n 



+ _ 



Hp<P 2 <p+l/L 2 ) + t{p 2 > p + l/L 2 ), 



and retain II . We get two terms Qq and . We estimate p^ in L°° as 
in (|5.10[) . getting 

IpM < ^ V#{p <k 2 <p + 1/L 2 } < p} < ^, 

since p < 1/L 2 (each of the two sets above contains a finite number of points which 
does not increase with L). Finally, we estimate p\ as in (|5.11[) and obtain 



\Tr LHCL) (VQt-)\< 
This time we have 



iru- 



iu 



- P 



1/2 



Tr 



L 2 (C L )\ 



-A L - p)Q- 



TTV-, 



Ai-Ml 1 / 3 



— y 



P 2 <M 
q 2 >p,+ l/L 2 



\V(p-q)\' 
q 2 - p 



<^ 2 E \V(k)\ 2 <CL 2 -* J |Vf 

q 2 >^+l/L 2 

Since d > 2, this completes the proof of Theorem 15. II 



□ 



5.2. Thermodynamic limit. With the Lieb-Thirring inequality (|5.2[) at hand, we 
can now relate the well-defined total relative energy in a large box to the one we 
have defined in Section 12.21 The following can therefore serve as an a posteriori 
justification of our definition of Try (—A + V — p)Qv- 
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Theorem 5.3 (Thermodynamic Limit, d > 2). We assume that d > 2 and \i > 0. 

Let V be a real-valued function in L 1 (R d ) n L 1+d / 2 (M. d ). Then we have 

(5.16) Try (-A + V - fi)Qv 
= ton (-Tr L 2 (CL) (-A L + V1 Cl - aO- + Tr i2(CL) (-A L - M )_ - M jf ^ , 

where the left side is defined in Definition \2. 6 A and — Al is the Laplacian with 
periodic boundary conditions on Cl — [~L/2, L/2) d . 



Sketch of the proof. We quickly explain the main steps of the proof, which proceeds 
by showing an upper and a lower bound. 

Let us fix a smooth finite-rank operator Q £ JC which we write in the form 

K ( 

Q = [lt/\ u i)( u j\+<lt,j l<K«il + ^/k)v%'l + 9i,j \ v i)( v j\ 

where {ui)f =l and (vi)fL 1 are orthonormal systems of the kernel and the range of 
respectively. The constraint that — IT < Q = Q* < 1 — LT is then only 
reflected in the coefficients </« (see [TU1 App. B] for an explicit representation 

of Q^j)- By assumption, the functions U{ and Uj are all smooth. Now we build 
from Q a test state Ql in the box Cl, by simply replacing the u^s and v,'s by 
orthonormal sequences (tt$ and (vi y L)fLi m i respectively, the kernel and the 

range of LT~ L . We can do this in such a way that 1c L Ui,L ~^ u ii ^C L v i,L — > Wi, 
Ic^Vu^i, — > Viti and Ic^V^x, — > V«j in L 2 (R d ), as Z, — > oo. One simple way to 
realize that is to periodize the functions as 

ke(2irZ/L) d 

and then to orthonormalize the so-obtained system. Similar arguments have already 
been used and detailed in [3]. The test state is then defined as 



K 

Ql ■= 



Yl (vt,?\ u i,L){ u 3,L\ + ( ii,j \ui,L){v hL \+ q i +\v l , L ){u^ L \ + q WJ \Vi, L ){Vj, L \j 



and it satisfies the constraint — H L < Ql < 1 — LT L by construction. We also 
have 



lim (Tr(-A L - n)Q L + / V PQl ) = Tr(— A ~n)Q+ f 
L ^°° V Jc L J Jr- 

Because we obviously have a variational principle in the box, 



VPQ. 



Al + Vt CL - /*)- + Tr L 2 (CL) (-A L - M )_ - M / 1/ 

inf Tr(-A L - /i)Q + / p Q y, 

-nr <Q<i-n7 .. Jn, 
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we deduce the upper bound 



limsup ( -Tr L 2 (Ci) (-A L + V1 Cl - (J.)- + Tr L 2 (Cl) (-A L - p)- - p / V 

L->oo \ JCl 

<Tr(-A-//)Q+ / V PQ . 

From the variational principle (|2.23|) in the whole space and the density of smooth 
finite-rank operators in /C, as stated in Lemma 13.21 we conclude that 

limsup ( -Tr L 2 (Ci) (-A L + V1 Cl - A*)- + Tr L 2 (Ct) (-A L - p)- - p \ V ) 

L^oo \ iCi / 

< Try(-A + ^-/i)Qy. 

In a second step we prove the reverse inequality, with the lim sup replaced by a 
liminf. We consider a sequence L n — > oo realizing this lim inf. Denoting by 

Q n = l(-A Ln + Vt CLn < At) - l(-A iB < m) 
the corresponding state, we know from our estimates that 



A Ln -^ 2 Q^\-A Ln 



/' 



1/2 



A, 



1 1/2 



©i(i 2 (Ci,„)) 

<C, 



< C, 



6 2 (L 2 (Ci,„)) 

and that lc i „A > Q„ = Pn + Pn where ||p^|| ioo < C//( d_1 )/ 2 (X n ) _1 and is bounded 
in L 2 (M d ) + L 1+2 / d (]R d ). Passing to weak limits, using V G L x (M d ) n L 1+d / 2 (R d ), 
we deduce that 

lhninf (-Tr LHCL ){-A L + Vl Cl - A*)- + Tr L 2 (CL) (-A L - M jf ^ 



> Tr(i?++ + JT 



?±± 



where l Cin p Q „ ^ p and t CLn \ - A Ln - p| 7 Q„| - A L „ - Mr /2 lc,„ R d 
Similarly, we have lc Ln Qn\ ~ A Ln - p\ x ^\ Chn S weakly-* in 6 2 (L 2 (R d )) and 
tc L Qn^-C L Q weakly-* in B. We now claim that Q £ K, pq = p, S = 
Q\- A -pi 1 "/ 2 , and i? ±± = | - A - p| 1/2 Q ±=t | - A - pi 1 / 2 . All this can be seen by 
testing against smooth functions of compact support, and we skip the details. We 
conclude that 



liminf -Tr L a, c )\ 



-A L + Vt CL - p)_ + Tr i2(Ci) (-Ai - p)_ 

> Tr v (-A + V - p)Q v 



by the variational principle (12.231) . This completes our sketch of the proof of The- 
orem EDS □ 
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5.3. Extension to positive temperature. In this section we extend our results 
to smooth partition functions, following [7 . This means we consider a smooth 
function / : R — > R tending to zero at infinity, and we look for a lower bound on 
the formal expression 

(5.17) Tr(f(-A + V)-f(-A)-f'(-A)V 

Our results above dealt with the function fo jfl (x) = — (x — Here we typically 
think of the free energy for a Fermi-Dirac distribution at positive temperature T 
and chemical potential /i, corresponding to 

(5.18) f T ^(x) = -Tlog (l + e~^' T ) , 

which converges to /o, M in the limit T 0. We will, however, be able to treat 
general functions /, provided they are concave and decay fast enough at infinity. 
The trick is to write / as an average of the reference functions /o,^ as 

(5.19) f(x)= [ (x-X).f'(X)dX, 

Jm. 

leading to the formal expression 
(5.20) 

Tr (/(-A + V)- /(-A) - f'(-A)v) = - jf Try (-A + V - A)Q A ,y /"(A) dX 
where 

Qx,v := 1(-A + V < A) - 1(-A < A). 
When / is concave, the integrand in the right side of (|5.20[) is > since Try (—A + 
V — X)Q\ v < 0, hence the integral always makes sense in R + U {+oo}. We may 
thus use this as a definition for the left side. In the following result we justify 
this formal calculation by a thermodynamic limit, and we state the corresponding 
Lieb-Thirring inequality. 

Theorem 5.4 (Lieb-Thirring inequality for smooth partition functions, d > 2). 
Let f : R — » R be a concave function such that f" € LJ^ C (M) and 

(5.21) / A 1+ ^|/"(A)|dA <oo 

Jo 

for some d>2. Then, for V £ L 1 (M. d ) n L 1+d / 2 (R d ), we have 

(5.22) L lim |Tr i2(Ci) f(-A L + V1 Cl ) - Tr L2(Ci) f(-A L ) - jf p n -* t) v} 

= - f Tr v (-A + V — X)Q\.y /"(A) dX, 

JR 

where, as before, — A^ is the Laplacian with periodic boundary conditions onCh = 
[—L/2,L/2) d . Moreover, we have the following inequality 

(5.23) /Tr y (-A + V- X)Q xy f"(X)dX 
Jm 

< L(d) J dXf"(X) J dx ((V(x) - A) 1 ^ - A+ + ^ + a| V(x)j . 
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The result holds the same under a weaker assumption on / than (|5.21[) , provided 
that the right side of (|5.22p is interpreted in a suitable manner. As such, Theo- 
rem [ITU already applies to the Fermi-Dirac free energy /t,^ as given in (|5.18p . since 
we have 

fruW = a 

' M r(i + e (^-M)/T) 2 

in this case. 

Proof. The Lieb-Thirring inequality (15.23)) is an immediate consequence of (I2.24[) 
and we only explain the thermodynamic limit (I5.22[) . First, it follows from the 
integral formula (I5.19[) and our assumption (|5.21[) . that /(— A^) and /(— + 
Vtc L ) are both trace-class. Using f|5 . 19[) . we obtain the identity 



Tr L2(C[) /( - A L + V1 Cl ) - Tr i2(Ci) f(-A L ) - f p ri -A L )V 

Jc L 

= - J^TT L 2 (CL) (~A L + Vtc L -X)Qx,v,l /"(A) dX 

with Qx,v,l = 1(-Ai + VI Cl < A) -t(-A L < A). By the Lieb-Thirring inequal- 
ity (|5.2p in the box, we have for A > 1/L 2 

- Tt lHCl) (-A l + V1 Cl - X)Qx,v,l 



d A 2 + d &„, , A 2 



< L(d) / [ (V(x) - \t +2 - A 1+ 2 + _±£ \2 V(x) + 2__\v(x)\ I dx 
Jc L 



( d-2 X 2 

< C 1 + A— + — — 



The last estimate is obtained by first replacing the domain of integration Cl by R 3 
(the integrand being > 0), and then using that 

((^A>^-A-f + ^n*)U ~ / \Vt 

\ 2 J A^OO 8 J R3 

For < A < 1/L 2 , we use (|5.3|) instead and obtain 

- Tt L 2 (Cl) (- A L + V1 Cl - X)Qx,v,l <C(1 + L- d ) . 

Finally, for A < 0, we simply note that 

- Tr L2(CL) (-A L + V1 Cl - X)Q x ,vx = Tr(-A L + V1 Cl - A)_. 

This last term vanishes when A < inf <r(— A^ + Vtc L ) an d ^ is bounded by 
Tr(— Al + V1c L )- < C(l + L~ d ) otherwise. As a conclusion, for L large enough 
we have a uniform bound 

(5.24) - Tr L 2 (Ct) (-A L + Vt Ch - A)Q W < C (l + A^) 1(A > —M) 

with M < limim C £_j. 00 inf cr(— A^ + Vtc L )- On the other hand we know by Theo- 
rem 15.31 that 

lim Tr L 2 ( c L )(-A L + V1 Cl - X)Qx,v,l = Tr v (-A + V - X)Q x ,v 

for every fixed A. Now (|5.22[) simply follows from Lebesgue's dominated convergence 
theorem. □ 
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6. Extension to more general background operators 

In the previous sections, we have considered perturbations of a constant density 
Pq. Our approach is, in fact, more general and we explain now how to handle 
other background densities. We typically think of a periodic background but, since 
we actually need very few assumptions, we state below an abstract theorem. Wc 
comment on the assumptions in the periodic case in Section [6.21 

6.1. An abstract Lieb-Thirring inequality with positive background. We 

consider a bounded-below self-adjoint operator H in L 2 (M. d ,C 9 ), with d > 2, and 
we fix a real number p 6 K. We assume that there is a constant C and an e > 
such that 

(Al) Pn\H-p\<E)(x) <C(E + E d/2 ^j for all E > and a.e. x £ M d ; 
(A2) 



1(M < H < p + e) t{n - e < H < p) 



|H-/i| 1/4 l^-Ai| 1/4 
(A3) p (x) := /t>i(H< M )(x) > e for a.e. x € ) 



<C\\V\\ L2 for all V e L 2 (R d ); 
6 2 



We define 



h-oo.jj,)^ ) anu n = 1 - it 



We emphasize that (Al) implies that P%ih=h) = 0: hence that p, is not an eigenvalue 
of H. With respect to the projections II~ and II+ we can decompose any bounded 
operator Q = Q ++ + Q~ + + + Q . Similarly as in Definition 12.11 we define 
the relative kinetic energy by 

Tr (H - p)Q = Tr\H- p\ 1/2 (Q ++ - Q—)\H - p\ 1/2 

for any bounded self-adjoint operator Q such that \H — ^,| 1/ ' 2 (3 ±± \H — p] 1 ^ 2 are 
trace-class. 

Theorem 6.1 (Abstract Lieb-Thirring inequality, density version, d > 2). We 
assume that the bounded-below self-adjoint operator H satisfies (Al)-(AS). LetQ be 

a self-adjoint operator such that — LT~ < Q < II + and such that \H — p\ Q ±± \ H — 

1 1/ 2 

p\ are trace-class. Then Q is locally trace-class and the corresponding density 
satisfies 

(6.1) p Q £ L 1+ 2(R d ) + L 2 (R d ). 

Moreover, there exists a positive constant K ( depending only on d, q, p, C and e ) 
such that 

(6.2) Tv (H-p)Q 

> K j ( (p (x) +p Q (x)) 1+d - Po {x) 1+ d - ——p { x )2p Q {x)\ dx 

with po(x) the background density ofTl~, defined above in (A3). 

Remark 6.1. For simplicity we restrict our attention to d > 2 but, with appro- 
priate modifications, a similar result holds for d = 1. If in Assumption (Al) the 
exponent d/2 is replaced by 5 > 1, our method still applies but the resulting lower 
bound is of course different. For instance, if relativistic effects are taken into ac- 
count, d/2 should be d in (Al), in which case the exponent 1 + 2/d in (|6.2p becomes 
1 + 1/d. 
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Remark 6.2. Similarly to Theorem I2.3[ it is possible to deduce from (|6.2|) a dual 
estimate on Tiy(H + V — p)Qv, where Qv = t(H + V < fi) — IT - , for any potential 
V G L 2 (R d ) n L 1+d/2 (R d ). For brevity we will not discuss this here. 

In applications, we typically think of H = —A + W(x) where W is a sufficiently 
regular function, and of p, strictly above the infimum of the essential spectrum of 
H. In Assumption (Al), the E d l 2 behavior of the density for large E is a rather 
general fact which we discuss below. On the other hand, the small E behavior in 
(Al) as well as (A2) are assumptions on H close to the Fermi surface. Vaguely 
speaking, (Al) is a (rather weak) assumption on the regularity of the spectral 
projections uniformly in x-space, whereas (A2) controls the interactions between 
particles inside and outside the Fermi sea. 

Next we show how to verify the large E behavior in (Al), under the assumption 
that H = —A + W(x) with W bounded from below. 

Lemma 6.1. Let W G Ll oc (R d ) such that W- G L°°(R d ), and consider the 
Friedrichs self-adjoint extension of —A + W on C^°(R d ). Then Pi(-a+w<e) * s 
uniformly bounded on M. d for every E £ M. and 

p [1(-A + W<E)}< (^) d/2 flW-Lco + E) d/2 , 
holds almost everywhere. 

Proof. Since W is bounded from below, we have, by the Feynman-Kac formula, 

r i pt\\W-\\ Lao 

P l \ - (4nt) d / 2 

where p[A](x) = A(x, x) denotes the density of an operator A. Using the inequality 
t(x < E) < e - t{ - x - E ^ we deduce that 

P [1(-A + W < E)] < e-p [e"*(-^] < { ^ )d/2 . 

Optimizing this bound with respect to t gives the result. □ 

6.2. Application to periodic backgrounds. In this section we restrict ourselves 
to periodic systems, that is, we take 

H = -A + W{x) 

where W is a Z d -periodic function which we assume to be sufficiently regular. Of 
course, we could as well consider other lattices than Z d . It is well known, see, 
e.g., [34j Sec. XIII. 16], that the spectrum of H is the union of bands 

a(H) = <j css (H) = |J |A„(£), £ G l-n,n} d } , 

n>l 

where A„(£) denotes the sequence of Bloch-Floquet eigenvalues of H with corre- 
sponding eigenvectors u n (£> %)■ Each X n is a periodic Lipschitz function of £, but 
the map £ u n (£) £ L 2 ((0,l) d ) is only piecewise smooth because of possible 
degeneracies. Writing for instance H = —A/2 + (—A/2 + W) > —A/2 — C and 
comparing the A n (£) with the eigenvalues of the periodic Laplacian in each Bloch 
sector, it can be seen that 

A„(0 >an 2 ' d -b 
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for some constants a, b > independent of £. Hence for every fixed /igE, there is 
only a finite number of n's such that A n (£) = fi for some £ S [— ir, ir) d . 
Let us fix fi > inf <7 eS s(— A + W). Then we have 



P o(x) = (27r)- d Y / I ^1(A„ 

„ >1 -/[-7r,7r)<* 



(£)<M)K(^)| 5 



Since ui(0, a;) is strictly positive, we easily conclude, by continuity in £, that po{x) > 
e > 0, and hence that (A3) is verified. 

Now we give some ideas on how one can verify Assumptions (Al) and (A2) in 
practice. First, we have 

PH\-A+W-„\<E)(x) = (2n)- d J2 [ d£l(|A„(0-H < E)\u n ^,x)\ 2 . 

Under suitable assumptions on W, u„(-,£) is bounded in L°°(M. d ), uniformly with 
respect to £, for each fixed n > 1. In this case, Assumption (A2) follows if the 
eigenvalues satisfy the following property: 

(6.3) | € [-7T, n] d , \X n (0 -n\<E}\< CE. 

This is generically true: If there is a unique n such that the graph of A„ crosses 
fi, and if V^A„(^) ^ for all £ with A„(£) = /i, one can easily verify that (I6.3[) 
is satisfied (/ii in Figure [TJ. At a point £ such that VjA„(^) = 0, the validity 
of (|6.3[) depends on the order of vanishing at this point. If, for instance, the second 
derivative is invertible, then ()6.3[) holds in any dimension d>2 (/i^ in Fig.QJ. Only 
the £'s which have a high (depending on the dimension d) order of vanishing can 
make (|6.3[) fail. When the Fermi surface is disconnected, each component being 
as before, the result is the same (^3 in Fig. [TJ. Finally, if A„(£) = A m (£) = /1 for 
n/m, the analysis is similar. For instance, transversal crossing of surfaces (/i4 in 
Fig. [TJ as well as Dirac-type cone singularities (/12 in Fig. [TJ are allowed. 

Verifying (A2) is much more subtle and requires a detailed analysis of the bands 
close to the Fermi surface. An exception is when fi lies in or at the edge of a gap, 
in which case (A2) is trivially satisfied (the estimate on pq±t in L 2 was already 
obtained in this case in [4]). In the case where /j, lies in the interior of the essential 
spectrum, we expect (A2) to be true, as soon as the Fermi surface is sufficiently 
regular. To make this intuition precise, a possible line of attack could be as follows. 
We assume again, for simplicity, that there is a unique n such that the graph of X n 
crosses \i. Then we have to prove that the operator whose kernel is 



-e<A n (f)<^ " Jn<Xn(Z')<H+e \j V ~ K {£,) \J K W) - l l 

is bounded on L 2 (R d ). The main idea is now that, for the question of boundedness, 
each Bloch function ti„(£, x) can be replaced by the corresponding plane wave 
exp(ia; • £). Arguments of this sort have been carried out in a similar context 
in [371 [3J[1 [S], for instance. When V 5 A n (£) ^ for all £ such that A„(£) = \i, the 
Fermi surface is smooth and can be locally replaced by a sphere. This reduces the 
computation to what we have done in Lemma [3T4l in the translation- invariant case. 

This concludes our intuitive description of how to verify Assumptions (Al) and 
(A2) for periodic backgrounds. Rendering all this rigorous is beyond the scope of 
this paper, however. 
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Figure 1. Typical Fermi surfaces for a two-dimensional periodic 
Schrodinger operator. 



6.3. Sketch of the proof of Theorem 16.11 The proof follows again the same 
two steps as that of Theorem 12.11 

Step 1. Estimate on Q ±± . We start by estimating the diagonal densities Pq±±- 
Following the proof of Lemma 13.31 Assumption (Al) implies that 



Tr\H- /i|7 > / R(p-y(x)) dx 



for every < 7 < 1, where 
R(p) = 

One has 

R{p) 



6Cq 



y/p- JCq{e + e d / 2 )) de. 



and R(p) 



00 (d + 4)(d + 2)(Cq) 2 / d 



P 



1+2 1 d 



Since e < p$ < M Yyy Assumptions (Al) and (A3), we therefore also have a lower 
bound 

R(p^(x)) dx 

- c J d ((Pv( x ) + Pi(. x )) 1+d ~ Po(x) 1+1 ~ ^-^P()(x)d p J (x)J dx 

with a constant c only depending on d, s, M and Cq. Applying these estimates to 
7 = ±(2 ±=t , we obtain the estimate analogous to (|3.9[) . 
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Step 2. Estimate on Q ±=F . Following the corresponding step in the proof of Theo- 
rem EH a bound of the form 



|p Q ± T | 2 <cTt q (H-^)Q 

<C\\V\\ L2 for &MV E L 2 (R d ) 



can be derived from the estimate 

n+ „ n- 



(A2') 



V 



|ff-Ml 1/4 I#-mI 1/4 

We now explain how to derive this bound from (Al) and (A2). 
To this end, we split 

IT = IJ- + HZ, = 10 - £ < H < n) + t(H < fj, - e), 

n+ = n+ + n+ = % < h < fi + s) + t(H > n + e) 

and estimate each term separately. The bound on II < VTI< is exactly our assump- 
tion (A2). For IL+VTI-, we write that 
(6.5) 



n: 



\H 



V- 



n- 



V4 \H 



1/4 



< 



1 



V- 



H 



1/4 



n- 



A 4 



1/2 



The density appearing on the right side is uniformly bounded on 
has, more generally, 

~t(\H -fi\ < a)' 



Indeed, one 



(6 ' 6) " \H-tf/* 

for every fixed a > 0. To see this, we write 
t(\x\ < a) _ 1 
jarj 1 / 2 ~ 2 
which, by (Al), implies the uniform bound 

"l(|if-Ml < a ) 



€ L° 



l(|a;| < min(y,o)) 



dy 

v 3/2 



\H~ / j,\ 1 /' 2 



C 
< — 
~ 2 



min(y, a) + min(?/, a) d ^ 2 



dy 

,,3/2 ■ 



Inserting (|6.6p in (|6.5[) and using the fact that H is bounded from below, we obtain 
the desired estimate for Tl^t V^II . The term corresponding to IT<Vn> is estimated 
similarly, using the fact that II> \H — /i] -1 / 4 < e^ 1 ^ 4 . 

This completes our sketch of the proof of Theorem 16.11 □ 
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